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1.  Introduction 

The  purpose  of  this  report  is  to  document  a set  of  computer  programs  for 
plasma  effects  on  one-dimensional  slot  antennas.  This  set  of  programs  covers  the 
cases  in  which  the  slot  is  located  in  an  infinite  conducting  plane,  an  infinite 
cylinder,  and  a sphere.  In  all  three  the  plasma  is  stratified  parallel  to  the 
given  surface.  Such  canonical  problems  are  approximate  but  useful  representations 
of  the  true  problem  of  plasma  effects  on  antennas  in  reentry  vehicles  and  interceptors, 
and  they  have  the  advantage  of  solvability.  In  applications  the  planar  and  cylindrical 
models  can  be  used  in  terms  of  the  tangent  to  the  vehicle  at  the  antenna  station 
and  the  spherical  model  is  useful  when  the  antenna  is  located  on  a spherical  aft 
dome  of  the  vehicle. 

The  types  of  plasma  effects  covered  by  these  models  include  the  antenna  radiation 
pattern,  the  input  impedance  and  the  antenna  noise  temperature.  The  models  are  a 
set  of  four  computer  programs.  SLOP  is  the  model  of  a slot  in  a ground  plane,  SLOC 
is  the  model  of  a slot  in  a cylinder,  SLOS  is  the  model  of  a slot  in  a sphere,  and 
ABCD  is  the  model  of  a two-port  network  which  is  used  to  derive  input  impedances  from 
the  results  of  the  former  programs.  All  are  written  in  FORTRAN  and  stored  on  the 
author's  user  file  in  the  GE,  Space  Division,  Information  Systems  and  Computer 
Center's  L66  computer. 

Section  2 gives  a general  description  of  the  models  in  terms  of  the  theoretical 
assumptions  and  techniques.  Section  3 gives  the  information  for  the  user  beginning 
with  input  and  output.  This  section  also  presents  a discussion  of  the  detailed 
theory  integrated  closely  with  the  FORTRAN  algorithms.  Generalizations  of  the  models 
which  are  possible  are  discussed  in  section  4.  Such  a generalization,  covering 
the  important  area  of  two-dimensional  aperture  antennas,  already  exists  in  a 
GS-RESD  computer  program  which  has  never  been  formally  documented.  The  latter  program 
is  quite  old  and  has  many  features  which  make  it  undesirable  for  future  revisions. 

The  present  models  are  thus  planned  for  such  future  revisions. 


2.  General  Description  of  the  Models 

These  models  assume  that  the  antenna  can  be  represented  physically  as  a one- 
dimensional slot  in  a perfectly  conducting  surface  which  is  either  a plane,  a 
cylinder,  or  a sphere.  The  plasma  is  represented  physically  as  a medium  external 
to  the  antenna  surface  having  its  properties  dependent  on  the  normal  distance  only. 
The  internal  representation  of  the  antenna  uses  the  assumption  of  a two-  port  linear 
network  where  one  port  is  at  the  input  terminals  and  the  other  is  the  slot. 


The  assumed  physical  representation  of  the  problem  makes  for  a convenient 
mathematical  solution.  The  dominant  mode  of  excitation  of  a one-dimensional  slot 
in  a conducting  surface  is  that  having  a uniform  electric  vector  across  the 
slot.  This  uniform  field  mode  is  assumed  to  exist  under  all  conditions.  In  the 
regions  external  to  the  antenna  surface  the  electromagnetic  fields  are  the  super- 
position of  separated  solutions  of  Maxwell's  equations.  For  the  planar  and  the 
cylindrical  geometries,  the  separation  constant  can  take  on  any  real  value,  so  the 
external  fields  are  represented  by  a Fourier  integral.  For  the  spherical  geometry 
the  separation  constant  can  have  only  the  discrete  values  associated  with  the  Legendre 
functions,  or  spherical  harmonics.  The  physical  meanings  of  these  separation 
constants  are  such  that  for  a given  value  the  em  fields  are  those  for  a plane  wave 
in  the  planar  case,  a cylindrical  wave  in  the  cylindrical  case,  and  a spherical  wave 
the  spherical  case.  The  effect  of  the  plasma  on  each  kind  of  wave  is  calculated 
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by  breaking  the  plasma  profile  into  a series  of  uniform  parallel  layers  and 
matching  boundary  conditions  at  each  layer  interface.  The  boundary  condition  at 
the  outer  surface  of  the  plasma  is  such  that  the  wave  is  outgoing,  or  radiating. 

At  the  antenna  surface  the  superposition  of  all  separated  waves  must  match  the 
assumed  field  distribution  in  the  slot. 

The  antenna  radiation  pattern,  including  the  plasma  effect,  is  calculated 
by  carrying  the  em  field  representations  described  above  into  the  asymptotic 
far  field.  For  the  planar  and  cylindrical  cases,  the  Fourier  integral  is  evalu- 
ated by  the  method  of  stationary  phase,  which  gives  the  result  in  terms  of  only 
the  wave  which  propagates  along  the  line  of  sight.  For  the  spherical  case  the 
far  field  involves  the  sum  of  spherical  waves  of  all  orders. 

The  near-field  effect  can  be  represented  in  terms  of  the  aperture  admittance, 
which  is  defined  as  the  ratio  of  the  tangential  magnetic  field  to  the  tangential 
electric  field  in  the  slot.  This  parameter  is  calculated  by  invoking  conservation 
of  energy  between  the  field  representations  in  the  aperture  and  in  the  medium 
outside  the  antenna  surface.  For  the  planar  and  cylindrical  cases  the  aperture 
admittance  is  the  integral  of  the  admittances  of  the  separated  waves  weighted  by 
the  square  of  the  Fourier  transform  of  the  aperture  electric  field.  For  the 
spherical  case  it  is  the  sum  of  the  admittances  of  the  separated  waves  weighted 
by  the  square  of  the  spherical  harmonics  expansion  coefficients  for  the  aperture 
electric  field. 

Noise  power  generated  by  the  plasma  and  received  by  the  antenna  is  calculated 
by  invoking  Kirchoff's  law.  This  law  is  used  in  terms  of  the  statement  that  the 
noise  power  at  the  antenna  input  contributed  by  an  element  of  volume  of  the  plasma 
is  equal  to  the  power  absorbed  by  that  element  of  volume  when  the  power  input  to  the 
antenna  is  equal  to  the  power  radiated  by  a black  body  at  the  temperature  of  the 
given  element  of  volume.  The  total  noise  power  is  the  sum  over  all  elements  of 
plasma  volume  of  such  noise  power  contributions.  Since  the  plasma  is  in  the  antenna 
near  field,  this  theory  gives  a result  which  is  similar  in  form  to  that  for  the 
aperture  admittance.  Thus  the  noise  temperature  at  the  aperture  is  a weighted 
integral  or  sum  of  the  noise  temperatures  for  the  separated  wave  solutions,  each 
of  which  is  a sum  of  the  noise  temperature  contributions  of  all  plasma  layers. 

The  transformation  from  the  aperture  surface  to  the  antenna  input  terminals 
is  done  by  the  program  ABCD,  which  utilizes  the  two-port  black  box  assumption.  It 
gives  the  relative  value  of  the  aperture  electric  field,  to  which  the  radiation 
pattern  calculation  is  normalized  in  the  other  programs,  both  with  and  without 
plasma  over  the  antenna.  It  also  gives  the  input  impedance  and  the  noise  temper- 
ature at  the  input  terminals. 


3.  Detailed  Information  for  the  User 

Section  3.1  contains  all  information  needed  to  operate  the  computer  programs. 
The  derivations  of  the  equations  are  Integrated  with  explanations  of  the  FORTRAN 
listings  in  section  3.2. 
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Figure  1 


Example  of  a Block  Data  Subroutine 


3.1  Input  and  Output 


Programs  SLOP,  SLOC  and  SLOS  have  similar  input  and  output  requirements,  rhree 
means  of  providing  input  are  used:  a block  data  subroutine,  NAMELIST  and  READ 
statements.  The  block  data  subroutine,  which  must  be  added  to  each  program  before 
it  will  run,  has  the  same  requirements  for  all  three.  An  example  is  given  in  fig.  1. 
This  subroutine  provides  the  data  on  the  plasma  profile.  COLL  is  the  collision 
frequency  in  units  of  10^/sec,  EMO  is  the  electron  density  in  cm“^»  TO  is  the 
temperature  in  °R  and  YO  is  the  distance  normal  to  the  wall.  NPTS  is  the  number  of 
points  in  the  profile  and  UNCM  is  the  unit  of  distance  in  cm  used  in  the  profile 
variable  YO.  For  example  if  the  data  YO  were  in  terms  of  Y/Rn  an&  RN  were  0.05 
inch,  UNCM  would  be  equal  to  0.127.  In  the  example  shown  the  profile  is  one  cm 
thick.  The  collision  frequency,  electron  density,  and  temperature  have  maximum  values 
of  unity  for  convenience  in  this  example  since,  as  we  will  see,  NAMELIST  input  permits 
scale  factors  to  be  used. 


1 


The  NAMELIST  name  is  INPUT  for  all  four  programs.  SLOP,  SLOC  and  SLOS  take 
basically  the  same  list  of  variables.  CAYA  is  k times  the  aperture  width,  which 
is  measured  along  the  surface  of  the  sphere  for  SLOS,  where  k is  2 7T  divided  by  the 
wavelength.  CYL  is  k times  the  body  radius.  In  the  planar  case  CYL  serves  the 
purpose  of  defining  a finite  length  for  the  aperture  so  the  output  noise  temper- 
ature will  be  finite.  In  this  case  the  aperture  length  is  27T  times  CYL,  as  though 
the  aperture  were  wrapped  around  a cylinder  of  radius  CYL.  Of  course  CYL  represents 
the  cylinder  radius  in  SLOC  and  the  sphere  radius  in  SLOS.  FACC,  FACE,  FACT  and 
FACY  are  scale  factors  for  the  profile  variables  collision  frequency,  electron 
density,  temperature  and  normal  distance,  respectively.  All  corresponding  values 
from  the  block  data  subroutine  are  multiplied  by  the  input  values  of  these  scale 
factors.  FMHZ  is  the  frequency  in  MHz.  GA  is  the  normalized  aperture  conductance, 
which  is  used  in  SLOP  and  SLOC  primarily  as  a control  on  whether  to  calculate  the 
aperture  admittance.  In  SLOS  only,  MODES  is  the  number  of  spherical  wave  modes 
to  be  used  and  THETO  is  the  angular  position  of  the  center  of  the  slot.  In  SLOP 
and  SLOC,  only  CAYA,  CYL  and  FMHZ  are  mandatory  inputs,  all  of  which  must  be 
greater  than  zero.  In  SLOS,  THETO  is  also  mandatory  and  greater  than  zero.  The 
programs  assume  that  the  scale  factors  are  unity  unless  input  otherwise.  Input 
nonzero  GA  in  SLOP  and  SLOC  only  if  the  aperture  admittance  calculation  is  not 
desired.  In  SLOS,  a value  of  MODES  (currently  20)  which  is  consistent  with  certain 
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dimension  statements  is  assumed  unless  input  smaller.  Do  not  input  MODES  larger 
than  this  value. 

Fig.  2 Shows  a sample  run  of  SLOP.  This  run  incorporates  the  block  data  in 
Fig.  1.  which  represent  a parabolic  electron  density  profile  which  goes  to  zero 
at  points  one  cm  apart  and  has  uniform  collision  frequency  and  temperature.  The 
first  input  during  execution  is  via  NAMELIST,  which  sets  values  for  the  aperture 
width,  the  body  radius,  the  profile  scale  factors  and  the  frequency.  The  first 
line  of  output  gives  the  result  of  the  first  iteration  in  the  calculation  of 
aperture  admittance  and  plasma  noise  temperature.  Reading  across,  the  first 
number  is  the  number  of  iterations,  the  next  two  are  the  real  and  imaginary  parts 
of  the  aperture  admittance  in  free  space,  the  next  pair  is  the  aperture  admittance 
in  plasma,  and  the  last  is  the  plasma  noise  temperature  in  °K  without  accounting 
for  aperture  area  or  conductance.  (As  noted  below  the  conventions  for  our  theo- 
retical derivations  are  such  that  both  complex  parts  of  the  aperture  admittance 
are  normalized  by  the  characteristic  admittance  of  free  space  and  the  standard 
electrical  engineering  convention  is  the  complex  conjugate  of  ours.  Thus,  for 
example,  the  susceptance  in  mho  is  the  negative  of  the  second  number  of  the  pair 
divided  by  376.7.)  Immediately  after  each  line  of  this  kind  of  output  there 
appears  another  where  the  program  reads  an  integer.  This  integer  must  have 

the  value  unity  if  further  iterations  of  the  admittance  and  noise  temperature 
calculations  are  to  be  done.  After  the  third  iteration  in  this  example,  this 
integer  is  input  equal  to  negative  unity,  which  stops  the  iterations.  (Zero  is 
another  option,  which  will  continue  the  iterations  but  only  with  respect  to  the 
aperture  admittance  in  free  space.)  Then  the  final  values  of  aperture  admittance 
and  noise  temperature  are  printed  on  a line  which  does  not  contain  the  iteration 
number.  The  aperture  susceptance  in  free  space  is  corrected  for  truncation  of 
the  iterations  and  the  noise  temperature  is  made  to  Include  the  effects  of  the 
aperture  area  and  conductance  on  this  line  of  output.  This  line  marks  a change  to 
a second  mode  of  operation. 

In  the  second  mode  of  operation  the  program  reads  a floating  point  number  at 
the  equal  sign.  This  input  denotes  the  far  field  line  of  sight,  or  incidence, 
angle  to  the  normal  in  degrees.  Thus  in  the  example  the  first  input  is  for  normal 
incidence.  The  next  line  of  output  gives  the  far  field  gains  and  losses  for  a 
constant  value  of  aperture  output  power.  The  first  number  is  the  attenuation  in 
dB,  where  a positive  number  indicates  a loss  of  signal  in  the  given  direction. 

The  second  number  is  the  gain  in  free  space  in  dB  and  the  last  number  is  the  gain 
in  plasma  in  dB.  (These  gains  and  losses  are  defined  in  greater  detail  below.) 
The  example  illustrates  the  method  of  generating  the  plasma  effect  on  the  radiation 
pattern.  If  the  input  incidence  angle  is  90,  or  greater,  the  program  reverts  to 
the  first  mode  of  input,  NAMELIST.  In  this  example  the  same  case  is  rerun  except 
that  the  peak  electron  density  is  1012  instead  of  1.24  x 10^.  Also,  input  of  GA 
not  zero  causes  the  program  to  skip  the  aperture  admittance  and  noise  temperature 
calculation.  The  next  time,  GA  Bust  again  be  input  since  the  program  sets  it 
equal  to  zero  Just  before  reading  NAMELIST.  Finally,  inputting  the  incidence  angle 
less  than  zero  causes  the  program  to  stop  executing. 

The  sample  run  shown  in  Fig.  2 illustrates  an  interesting  phenomenon  in  non- 
uniform,  relatively  collisionless  plasma.  That  is,  the  attenuation  at  large 
incidence  angles  for  the  electric  vector  in  the  plane  of  incidence  is  maximum 
when  the  peak  electron  density  is  near  critical  density.  (Critical  density  in 
cm  -3  is  approximately  equal  to  1.24  x 10-8  f2f  where  f is  the  frequency  in  Hz.) 
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Fig.  3 shows  a sample  run  of  SLOC.  All  Input  and  output  formats  are  identical 
with  those  for  SLOP.  Even  the  numerical  results  are  qualitatively,  although  of 
course  non  quantitatively,  similar. 

Fig.  4 shows  a sample  run  of  SLOS  for  the  same  case,  but  where  the  slot  is 
centered  at  45°  from  the  pole  of  the  sphere.  In  the  iterated  output  of  aperture 
admittance  and  noise  temperature  the  integer  denotes  the  number  of  modes  used. 

The  iteration  is  continued,  up  to  a maximum  of  20  times  by  inputting  unity.  It 
is  stopped  by  inputting  zero.  There  are  no  other  alternatives  since  it  is  logically 
difficult  to  divorce  the  calculations  with  and  without  plasma  in  the  spherical 
case.  The  ouput  of  the  final  results  shows  the  total  number  of  modes  used  because 
there  is  an  internal  criterion  on  truncation  of  the  number  of  modes.  The  operator 
must  watch  to  see  that  the  number  of  modes  continues  to  increase  when  he  inputs 
unity.  If  the  output  number  should  be  the  same  as  on  the  previous  output  line 
or  if  the  number  20  is  reached,  the  program  has  switched  from  this  mode  of  input/ 
output.  Then  the  input  consists  of  the  polar  angle  of  incidence  in  degrees  for 
output  of  the  far  field  attenuation  and  gains  without  or  with  plasma.  For  this 
sample  case  the  normal  incidence  direction  is  45° f the  grazing  incidence  direction 
is  135°  and  179°  is  deep  within  the  optical  shadow.  Input  of  180°,  which  is  a 
singular  point  for  the  assumed  antenna,  switches  the  program  back  to  the  original 
input/output  mode.  Input  MODES  less  than  20  to  skip  the  iteration  of  the  output 
of  admittances  and  noise  temperature.  Then  the  program  goes  immediately  to  the 
far  field  pattern  mode  of  input/output.  As  shown  by  this  example,  MODES  does 
not  have  to  be  input  again  to  repeat  this  modus  operand! . To  stop  execution 
input  0.  when  the  program  reads  the  incidence  angle.  Note  that  there  is  no  option 
to  skip  the  part  of  the  program  which  calculates  admittances  and  noise  temperature, 
because  this  is  a logical  impossibility  in  the  spherical  geometry. 

Fig.  5 illustrates  a run  of  the  program  ABCD.  The  NAMELIST  inputs  are  the 
variables  P,  R,  Y,  YAO,  YO  and  YY.  P and  R are  arrays  of  phase  (in  degrees) 
and  amplitude,  respectively,  of  the  (experimentally  measured)  complex  voltage 
reflection  coefficient  at  the  antenna  input  terminals.  For  each  of  these  the 
first  element  is  for  the  antenna  in  the  free  space  environment,  the  second  is  for 
a good  conductor  tightly  covering  the  antenna  aperature,  and  the  third  is  for 
a thin  resistance  sheet  covering  the  aperture.  The  complex  array  Y is  the  input 
admittance  corresponding  to  each  of  these  same  three  conditions  in  the  same  order. 

Of  course  Y is  redundant  with  P and  R,  and  it  represents  an  alternative  method  of 
input.  If  the  phase  and  amplitude  of  the  reflection  coefficient  at  the  antenna 
input  are  input  to  the  calculation  then  Y need  not  be  input.  If  it  is  desired  to 
input  Y rather  than  P and  R,  P(3)  must  be  input  equal  to  a negative  number,  or 
set  P-3*-l.  for  example. 
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YAO  Is  the  complex  aperture  admittance  In  free  space.  YO  Is  the  characteristic 
admittance  of  the  Input  line  and  YY  Is  the  conductance  of  the  thin  resistance 
sheet  which  corresponds  to  the  third  set  of  data  for  the  complex  reflection 
coefficient  or  admittance  at  the  input  port.  The  units  of  Y,  YAO,  YO  and  YY 
must  be  consistent  and  the  phase  convention  for  P must  be  consistent  with  the  complex 
variable  convention  for  YAO.  Since  a lossless  input  line  and  an  ideal  resistance 
sheet  are  assumed,  the  variables  YO  and  YY  are  real.  By  default,  as  in  the  example 
shown,  YO  corresponds  to  a 50  ohm  line,  YY  corresponds  to  a 100  ohm  resistance 
sheet,  and  the  characteristic  admittance  of  free  space  is  the  unit  of  admittance. 
Except  for  the  values  of  P and  R,  which  have  been  assumed  arbritrarily , the  data 
for  this  run  have  come  from  the  results  of  the  runs  from  figs  2 through  4 in  order. 

The  first  eight  numbers  output  are  the  complex  values  of  the  A,  B,  C and  D 
parameters,  which  are  defined  in  section  3.2.4.  The  next  output,  on  a line  by 
itself,  is  the  reflection  loss,  1-|r|2,  in  dB  for  the  antenna  in  free  space. 

Then  the  program  reads  input  values  of  the  complex  aperture  admittance  in  plasma 
and  the  plasma  noise  temperature  at  the  aperture.  Qf course  the  units  for  the 
aperture  admittance  must  be  consistent  with  those  for  YAO.  The  first  two  numbers 
in  the  next  line  of  output  are  the  amplitude  and  phase  of  the  voltage  reflection 
coefficient  at  the  input  when  plasma  covers  the  antenna.  Next  is  the  correspond- 
ing reflection  loss  in  dB.  The  fourth  output  is  the  angle-independent  part  of  the 
total  radiated  signal  loss  in  dB.  This  latter  number  must  be  added  to  the  attenuation 
numbers  output  by  the  other  programs  as  a function  of  angle  in  order  to  give  the 
total  signal  attenuation  relative  to  the  signal  in  free  space,  which  includes 
reflection,  absorption  and  antenna  pattern  distortion  effects.  The  last  output  is 
the  noise  temperature  at  the  antenna  input.  As  shown,  one  normally  iterates  input 
values  of  admittance  and  noise  temperature  at  the  aperture,  which  allows  the 
operator  to  repeat  the  calculations  conveniently  for  a number  of  different  plasma 
conditions.  Inputting  zero  noise  temperature  causes  the  program  to  go  back  to  new 
NAMELIST  input  and  inputting  negative  noise  temperature  makes  it  stop  executing. 

The  third  set  of  NAMELIST  input  illustrates  input  of  Y instead  of  P and  R. 


3.2  Computational  Techniques 

This  section  of  the  report  gives  the  theoretical  derivations  as  well  as  the 
explanation  of  the  FORTRAN.  These  derivations  use  a peculiar  system  of  units, 
which  is  designed  for  maximum  convenience.  The  unit  of  length  is  the  free-space 
wavelength  divided  by  2TT.  The  units  of  the  electric  field  E and  the  magnetic 
field  H are  the  same;  l.e.,  the  characteristic  impedance  of  free  space  is  unity  and 
all  admittances  in  our  equations  are  numerlcallyequal  to  their  values  in  mho 
times  376.7.  All  fields  have  the  time  variation  exp  (-iajt) , which  is  suppressed 
in  all  equations.  Thus  a wave  traveling  in  the  positive  x direction  has  the  vari- 
ation exp  (ix).  Another  consequence  of  this  convention  for  complex  quantities 
is  that  the  standard  electrical  engineering  results  are  equal  to  the  complex 
conjugate  of  these  results. 

The  programming  techniques  are  straightforward.  Most  of  the  built-in  functions 
and  library  routines  which  are  called  are  of  the  every  day  variety.  The  only 
exception  is  FXOPT,  which  is  used  only  to  suppress  the  error  message  for  exponent 
underflow.  (The  L66  computer  word  does  not  allow  numbers  smaller  than  about  10-38; 
and  when  such  a number  is  generated  during  execution,  the  computer  sets  it  to 
zero  and  normally  prints  an  error  message,  which  FXOPT  is  used  here  to  suppress.) 
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20  COKkOh/H  AC/'l  ( Ovj ) ,YA(4) 

Ju  COfr.PLh  A CAtCbtCC,L)IbL(Oo),U(4)  ,YA,  YAO 

40  UI  it.cNblOw  Ch  IT(  2 ) , dM(  2 1 ) ,X(  Ov)  , AA(  o)  , Y ( 2 I ) 

‘jo  hA/».t:LI  bl  /I  uPOT/LaY  A, Cy  L,  FACC,  f ACb « FACT , FACY , FMriZ.UA 

00  DA'iA  FACC, FACb, PACT, FACY/4*!  ./ 

70  uA  1 A P I , kAU/3.  1 4 1 b927 , . v,  I74b329/ 

bv  CALL  PXOPT(o/,  I , I ,v,) 

Vu  1 JA-o . 

1 \~i\j  ithAL  IbPui 

I 1 L=  I 24L-J . I b*F/.  hZ*FMllZ/FACb 

I2u  CoLA=.  I 59  1 b494>.  FACC/Fi/tliZ 

I Jo  CAY  =2 . v,9bb44  /b-4*  Ui4Ck*Pi/iHZ 

1 4 u Child)-). 

I bo.  Jit  1T(2)  = I. 

1 00  3 A=  I 

IVv,  JO  lu  I»lfMPTb 

lov,  Y ( I )=FACY*YO(  I ) 

I 9 v,  l o,  Fj-i  ( I ) =btoO  ( 1 ) /Fa  o 

200  CALL  Ch  IT b ( bM , \ , hPTb , CH  1 1 , XA , A ) 

2lv,  a=aA(  OA  )-aA(  1 ) 

2^v,  iJ=0 

^ Jo  mP^Z 

24o  JO  00  I=2,0A 

2b O J=i-I 

200  J-AA(J) 

2/0  J— A A ( I ) -U 

2oo  A -bO.<to/A 

2Vo>  I F ( u . L 1 . 1 ) h=  1 

Joo  0 ~ U /Y\ 

Jlv,  b-0AY*C 

32o  JO  bO  L=  I . K 

JOO  .4— *i+  I 

040  X(U)=b 

jbo  F=U+(L-.b )*C 

Jou  Jo*  Jv,  M—I»\P  , NPTb 

37„  IF(Y(i/i).Ol.P)0l,  TO  35 

3bo»  JO  Cu«!InUb 
39  0 Jb  MP=m 
400  WWi—lflP- 1 

410  AA=(F-Y(Mrt))/(\ (mP)-YUm)  ) 

4^o,  bb=  1 . -AA 

43v,  COL-COLA* (C0LL(  MM)*t>b+COLL(Mp)*AA) 

44  v,  i'(i4 ) =/•  ACi  *(T0(  h i<i  )*bd+ AA*TO(  MF  ) ) 

4bv,  IF(  AA.oc.  I . >00  TO  4o 

4oo  I F ( h*M ( tart)  • OT.  o.  . AND  .Hm(  mP  ) • Ui’  ,0 • ) 00  10  3b 

4/  \j  0=-b.<i(  i.ii.i)*bb+tM(  |.\P  )*AA 

400  00  10  4d 

49  v,  Jb  J=bi<t(  i.nvi  )**Bb*bi<.  ( mP  )**  A A 

bLO  00  10  4b 
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Figure  6 Listing  of  Program  SLOP 
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LI8T 

b 1 0- 1 oOv/ 

blu 

4<j 

J=-bi<i(  i>tH ) 

b2o 

4b 

J*0/(  1 . +C0L*C0L ) 

b3o 

a-o*00L 

b4  0 

bw 

GltL<U)=CmPLX( 1 .-0, H) 

bb  kj 

60 

Cun  YI  ilUt 

bo 

IbCOA.OI. 0. )00  TO  dO 

b7G 

0 — VJ  • 

bbO 

M=l 

b9U 

ITtK=Vy 

OOU 

CALL  HkOP(l>i,DltL,X,M,U,Q) 

Oi  vj 

1 1-( U.LT.o. )G0  70  6w 

02u 

AA=.9v,28 

Oo  v 

YAO=CmPLX(AA,u. ) 

6 4U 

C0=(0(  l)-u(3))/(0(4.)-Q(2)) 

6bw 

Ya=AA*CC 

OOu 

A=«hAL(CC) 

67  U 

O-Al^AOtCC ) 

obu 

Ti«=AA*(  TA(  1 )+A*  l"A(  2)+b*TA( 

o9G 

D0=2.  *P  I/CAYA 

7 uu 

U=1 .b*DU 

7lu 

■ik=  1 . o 

720 

L-  i 

7 ju 

61 

CALL  PkOP(W#l.)It  Lf  X,M,  U,G) 

lA\j 

IMU.Ll.v,. )L=0 

7 b(J 

o2 

ou=l  ./(Pl*(bA-.Ob/liA) )**2 

7 Ou 

i>0=  1 . -U*U 

770 

Ib(60)6b,69,o3 

7bO 

0 3 

a-bOHi(bC) 

79u 

YAO =Y  AO+do/h 

600 

Ib(L.«£.l)UJ  10  69 

610 

CC=(U(1  )-n*0( 3 ) ) /( >1*0  ( 4 )~0 

020 

GO  TO  66 

6 3<j 

6b 

»<=80hT(-6C) 

840 

YAO"YAO-CMPLX(v. , 6d/h ) 

6bO 

Ib(L.ilH.  1 )00  TO  69 

860 

CA=CfoPLX( u< t M ) 

8/0 

CC=(0( 1 )-CA*G ( 3 ) ) /( CA*0(4 ) 

6bo 

ofc 

i a=YA+OU*CC 

690 

A=KbAL( CC) 

900 

d=A  IMA0(  CC) 

9 1 u 

Tl«=Tw+Bb* ( i A(  1 / +A* 1 A(  2)  +0* 

920 

69 

8A=  6A+ 1 . 

93o 

0*0+00 

94  o 

ritk*ITLri+1 

9b0 

Phial  1 va^o  « 1 i bh  « YAO  * Y A « fu 

960 

hbADiL 

97  o 

I h( L) 70,62,6 1 

96u 

70 

A®. 2b*CAYA/PI /( PI*( OA-.b ) ) 

99w 

YaO=Y AO+CmPLa ( 0 . , -A ) 

1 woO 

OA*hhAL(YA) 

> ) + ( A* A+  b*b ) *TA ( 4 ) ) 


Figure  6 (continued) 
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lolL 

lN=iN*CAYA*CYL/GA 

1 okO 

FKIIVI*  3oOO,YA(.  ,YA,  fn 

CoO 

bo 

kbAlUThb'l 

lkTlnbi.Ll.O.)  00  10  100 

0 jO 

IMTHfcT.GE.9o.  )U0  TO  1 

1 uOu 

L=b  1 N ( HAD*  lHb’j  ) 

o /o 

CALL  PH()P(W.UlbL,XtL,0,U) 

OoU 

IbCO.LT.o.  )00  TO  do 

o 9o 

In=b0xl  ( 1 . -U*U) 

1 Oo 

L=.b*CAYA/UA . 

1 1 10 

0=  .b*CAYA*U 

1 1 4o 

lb  (D.UT. o . ) 0=0* ( o i N.(  0)  /J  ) **k 

1 jo 

0=4. J43*AL0O(C) 

1 1 4o 

CA=0(  4 )-Q(  k)  /\. 

1 1 bo 

Cb=0( 1 )*Q  ( 4)  -0  ( k ) *C(  3 1 

1 loo 

CC=CA/Cb 

1 1 70 

A=4 . j4  J*  ALOU  ( k EA  L ( COCON  JO  ( CC ) ) ) 

1 1 bo 

0=0- A 

1 1 Vo 

pH  INI  kooO , A , L , C 

1 k oo 

00  TO  bo 

Ik  lo 

1 00 

bTOP 

1 kko 

1 OO0 

kOKkAi  ( lo, jX,  Ikkb  Ik.  4,  jX  , 1 Pkb  1 2. 4 ,.|  PE  1 b.  4 ) 

l4j0 

koOU 

k0kkAT(bk9.k) 

Ip  4o 

300o 

kOhiv.Ai  (ox  , 1 Pkb  1 k . 4 , oX , 1 Pkb  1 k.  4 , 1 Pb  1 o . 4 ) 

Ik  bo 

tziiD 

Ikoo 

bObnoLi  ii«b  PhCP(  h , ul  bL , X , i i , V , 0 ) 

1 4 1 o 

CWi-.mON/ NAC/' 1 ( £ 0)  , 1 A ( 4 ) 

1 k bO 

COniPLbX  CA , CO , CC  t CO,  Cb  ,Cr  ,01  bH  oO ) , 3 ( 4 ) 

lk9o 

UlMbublOn  T Li  < ^ ) , lC(4  ) , A(  Oo) 

1 j oo 

\b=\ *V 

1 J 1 o 

Ik  (w.LT.o)Ol)  ’10  lo 

1 jko 

ll 

-3 

1 JJO 

T A ( I ) =0. 

1 34o 

b 

T b ( I )=o. 

I30o 

1 b ( k ) = 1 . 

1 job 

lo 

0(  l)=(  1 . ,o. ) 

1 37o 

0( k) = ( O. ,0.) 

1 jt)o 

0( 3) =Q (k  ) 

l39o 

0(4)=C( 1 ) 

1 400 

00  90  1=  1 ,N 

14  lo 

CA=CSGkT(DIEL(  I )-C..iPLX(  Vo,o.  )) 

1 4ko 

u=X ( I ) *A I MAO ( 0 A ) 

14  o0 

I k (Aob ( u) .CT. 34. b4)00  TO  loo 

14  40 

A=X(  1 )*NbAL(CX  ) 

14  bo 

o=.b*C03(A) 

|4  Oo 

u= .b*o IN ( A ) 

14  70 

b=bXP(b) 

1 4 bO 

k=l ./b 

l4  9o 

0=E+k 

lb  00 

h=l — b 

* 


Figure  6 (continued) 


I 


SHIS  PACE  IS  BEST  QUALITY  PRACTICABtl 

fBtfl  QQfcY  umausHHD  TODO0  - - -14- 


i_ibi  I b I o—«ioUo 


I L i v 

OL.-C.'  o' La  ( OJ , i.  «ii ) 

1 

1 b S.\u 

oo-L u'L.w  ( b*ri , two ) . 

i 

ib  bo 

'■#L* 

i b *>  o 

bo-uo/bA 

J 

1 *-  Jo 

bn— bu/Li fcL ( I ) 

1 J Ou 

br-bb^'Ll  t:L  ( I ) 

• 

1 -J  1 0 

bn-b.-‘*0(  1 )+GL;v  j) 

1 

1 1/  oo 

L ( o) -Co*b(  o)  +bT»Aj(  1 ) 

i 

- J 

lb  Vo 

L ( 1 ) -bn 

1 

1 t’Oo 

bn  — b b ) ^*b U ( 4 ) 

1c  lu 

b ( 4 ) =bt>fo(  •> ) -rbr*U(  ,:) 

lb  c'U 

L ( k ) =La 

1 

1 C JU 

i r ( i«  • L i • o ) uO  ib1  130 

\ 

1 1;  4o 

bi  — bOl i Jo  ( o( J ) ) 

1 OJt, 

bb=bb/i  Jo  ( u ( 4 ) ) 

\ 

1 ( Ob 

Toil  )-iitinL('j(  i )*olj) 

i 

l<-  To 

7 b ( k ) -niiALC  0 ( 1 ) *bC+G  ( k ) -fbl. ) 

J 

1 (.  Lb 

7b(  u) -Al.-iAuC  Gl  1 ) .Tb'C-u  l *Cl;J 

1 

1 O VO 

1 b ( 4 ) — iili/vL.  ( O ( £ ) A bb  ) 

i 

i / oo 

LO  no  J—  i , 4 

J 

1/  lu 

i A ( J ) =TA  ( J ) ■*•  l ( 1 ) * ( XL  ( j ) -7  b ( J ) ) 

1 1 I'b 

li.(  j)  = xb(  J) 

1 

i ■/  Ob 

'JW 

A— o . 

1 / AL 

I/O  »U  J—  1 , 4 

i 7 ‘Jo 

OU 

>.=A+ntAi-(  COoJb  ( Q ( J ) j *u(J  ) > 

1 I 

1 7 Oo 

Irtrt.ui.  I.tikioo  Yu  1 bo 

1 7 7u 

Vu 

Ob  l.  l I , i uu 

i '/  liu 

i.  r;  1 bUO 

« 

1 7 9o 

i 

\ — — » — 1 • K— bO 

1 t.Ut 

nhTi  n»; 

1 io 

i.i « ij 

i i bo 

b U f.n b L i I i i b wl£  1 Ib  ( 1"  I , a 1 , n , a « a A , nn  ) 

it'  bo 

l.  1 i.ihi/b i Ol i A(  k;  ,1-  i ( lO  t bi  ( o ) ,XA(  6)  , XI  ( U ) 

1 , 

1 C/4o 

On l A BTb/l.n— c/ 

1 1 

If.  Jo 

IT-iiA 

1 i:  ob 

i/A=  1 

■ 

1 o/b 

LO  1 1-1,4 

1 t oo 

i 

Xn  ( I ) — o . 

U Vo 

1 i ^ 1<  “*  i 

1 ‘ ub 

ub  lob  1-1,  IT 

IV  io 

Ob  9o  j — 1 « m 

iV<.'v> 

lv-  J+  i 

1 V Oo 

hi  .-l- 1 ( J ) /A ( 1 ) — 1 • 

IVAo 

Tr'=sl‘i  ( a)  /A(  I ) - 1. 

i V 'jo 

1: 

1 V 01. 

IT (tt.oi.U. )oO  TO  VU 

IV/o 

i »- ( t . ol.o.  )uo  To  j 

1 V 'JO 

1i  (ai  b(i-/..)  .J7 .0.  )uo  i L .hi 

1 Vvl 

a A (ha) “X  I ( J ) 

1 

('oJl- 

1 i*  ( !•  A • fcG , 4 )bo  TO  1 A- 

4 

Vf 

Figure  6 (continued) 

-15-  IHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

PROM  OQPX  JUKKISHED  TO  DOCS  - *" 

LlO'i  2uM-*:j6U 


20  1 v/ 

i,A=NA+  1 

oO  10  9xj 

c.\j  oO 

0 

Xi*i=XI  ( J) 

2040 

Xp=XI (K) 

tL\j  L)  w 

IE  ( Em.  Lh . - 1 . . LA . EH  . LE . -1  . ) 00  TO  10 

2060 

0=X-.+  (XH-XM)*/L0G(A(  D/FKJ)  )/ALOU(Ei(K)/HI( J)  ) 

20  / u 

00  TO  4u 

2uO(j 

10 

t=Xl>.+  ( XP— Xi<0  *EM/(  Em— HP  ) 

2v>9o 

40 

XA  ( i,  A ) =b 

2100 

li-(uA.EQ.4)GO  TO  120 

2110 

OA=OA+ 1 

2 1 2u 

90 

CONTINUE 

2 1 30 

1 00 

00 M I NuE 

21  40 

LO  110  1=2,  .Vi 

2 1 DO 

IE  (El  ( i-1  ) .LT.EI  ( I ) .AI.U.EKI ) .OE.EI  ( 1+ 1 ) )00  TO 

^ 1 Ow 

1 lo 

CONTINUE 

21  7o 

KA=UA-I 

1 Ovj 

00  10  I2u 

4 1 90 

1 1 0 

aA(imA)=XI(  I) 

*?2  uu 

120 

1(A=hA+  1 

22  10 

XA(hA)=XI(I) 

tLS  2w 

uA=i,A+  1 

222o 

Art ( LA ) =X  I ( ii) 

224o 

00  1 ju  1=1,  OA 

22  30 

b= 1 • u30 

,cT  Ou 

OO  1 4u  J — 1 , l,A 

2270 

iE(XA(J).OE.b)GO  TO  I4u 

2260 

i\—J 

<L/.  9 vj 

L=XA( J ) 

2300 

1 4 u 

COM  TnUE 

^ j 1 

01(1) =0 

2o  2vj 

1 00 

Art (K ) = 2. E30 

2o30 

00  1 60  1=1,1,  A 

234  0 

1 Oo 

XA ( I ) =01 ( I ) 

23 

hETllKH 

<ijOu 

Ei,U 

* 

Figure  6 (continued) 

(continued) 


-16- 


3.2.1  SLOP 

Fig.  6 Is  a listing  of  this  program,  which  consists  of  a main  program  and 
two  subroutines.  Of  course  the  block  data  subroutine  must  be  added,  as  discussed 
above . 


3. 2. 1.1  Main  Program 


After  reading  the  NAMELIST  type  of  input,  the  program  generates,  in  lines  110 
through  550,  a representation  of  the  plasma  profile  by  a series  of  uniform  layers 
in  terms  of  an  array,  DIEL,  of  the  complex  dielectric  constant  and  an  array,  X, 
of  the  layer  thickness.  The  equation  for  the  complex  dielectric  constant  K is 


K - 1 - (Ne/Nc)/(1+  i V/co) 


(1) 

where  Ne  is  the  electron  density  and  Nc  is  the  critical  electron  density.  Except 
for  the  electron  density  scale  factor,  FACE,  line  110  gives  the  equation 
for  the  critical  density.  The  array  EM  represents  Ne/Nc  at  input  plasma  profile 
points  and  the  array  Y represents  the  input  profile  distances  in  cm.  After  calling 
CRITS  (see  below),  there  are  NA  points  in  this  profile,  located  at  the  distances 
XA,  between  which  interpolated  uniform  layers  are  desirable.  The  algorithm  from 
line  210  through  550  performs  the  interpolation  accordingly,  with  the  additional 
constraints  that  the  total  number  of  layers  be  approximately  fifty  and  that  all 
their  thicknesses  be  as  nearly  uniform  as  possible.  The  collision  frequency  is 
linearly  interpolated,  as  is  the  temperature  array  T.  The  electron  density  is 
logarithmically  interpolated  where  possible. 


The  aperture  admittance  and  plasma  noise  temperature  are  calculated  in  lines 
570  through  1010.  To  derive  the  equations  for  these  quantities,  let  the  plane 
z-o  be  the  ground  plane  and  let  the  slot  be  between  the  lines  x = - a/2,  where 
a is  the  width  of  the  slot.  Then  in  this  plane  let  Ey  = 0,  let  Ex  = 0 if  |x|>a/2, 
and  let  Ex  * E0  if  |xKa/2.  Of  course  E0  is  the  aperture  field  strength,  the 
determination  of  which  can  be  left  for  later.  This  assumed  distribution  of  electric 
field  gives  rise  to  a magnetic  field  vector  which  is  everywhere  parallel  to  the  y 
axis,  so  H ■ L and  we  need  not  denote  a vector  component  as  witn  E.  Also,  both 
E and  H are  independent  of  y.  Maxwell's  equations  in  tins  situation  are: 


3 f*x  " 3E  z 
T-,  TT—  = iH 

a z 3 x 


(2) 


3 H 
3z 


iKEx 


(3) 


9 H = - 1KEZ  (4) 

9 x 


Accounting  for  the  fact  that  the  gradient  of  K is  parallel  to  z,  these  equations 
combine  to  give 

3 H + _ _i  ^K_  + kh  - 0 

K 3z  3z  (5) 

The  variables  separate  if  the  first  term  is  equal  to  a constant  times  H.  In 
order  to  give  solutions  which  are  finite  for  all  values  of  x,  let 


■ - u^h 

3x2 


(6) 


where  u is  any  real  number.  Since  the  medium  above  the  groud  plane  is  infinite 
and  since  there  is  not  an  infinite  number  of  periodically  located  slots  in 
the  ground  plane,  the  allowed  values  of  u are  continuous.  Therefore  the  principle 
of  superposition  requires  that  the  total  solution  be  an  integral  over  all 
solutions  of  (6);  i.e.. 


w 

H (x , z ) * J h (u , z ) exp  ( i u *)  du 


where  h(u,  z ) is  the  part  of  the  separated  solution  which  depends  on  z , and  of 
course  also  on  u.  This  equation  is  a Fourier  transform  and  its  inverse  is 


h (u,z)  - 7 Jr  /H<x»  z)  exP  (-  i u x)  dx 


A precisely  similar  pair  of  equations  can  be  written  for  each  component  of  E, 
so  the  vector  E(x,  z)  has  a transform  in  terms  of  the  vector  e (u,  z). 

The  aperture  admittance  Ya  is  the  ratio  of  the  magnetic  field  to  the 
tangential  electric  field  in  the  slot.  The  fact  that  this  quantity  is  treated 
as  a constant  implies  that  the  magnetic  field  in  the  slot  is  a constant,  as 
is  the  electric  field.  This  is  impossible,  so  let  us  use  the  conservation  of 
energy  to  estimate  Ya.  Thus,  using  * to  denote  complex  conjugate, 

oo 

a No  I 2ya  * J Ex  (x»  °)  H (x>  °)  dx  (9) 

where  both  sides  are  the  integral  of  E*  XH  over  the  ground  plane,  the  left- 
hand  side  using  H = Ya  EQ  in  the  slot.  Substitution  of  (7)  in  (9)  and  the 
assumption  that  the  order  of  integration  can  be  reversed  gives,  after  carrying 
out  the  integral  on  x as  the  transform  of  EXj 


a|Eo  | 2Ya  = 27r  J h (u,  o)  e*  (u ,o)  du 


Now  ex  (u,o)  is  the  transform  of  the  assumed  tangential  electric  field  dis- 
tribution in  the  ground  plane. 

d/2 

ex  (u,o)  = -|2p-  J"  exP  ( - i u x)  dx 

-a/z. 


J exp  ( - i u x)  dx 


E0  sin  (u  a/2)  / ( TT  u) 


Therefore 


Ya  “T 7* 


co  o 

2 f y(u)  [j'hfMi/aW f du 

rr  A J 


where  A u)  is  the  admittance  in  the  plane  z = 0 for  a wave  having  the  separation 
constant  -u  , or  having  exp  (i  u x)  variation  with  x. 

A wave  havingexp  (i  ux)  variation  has  exp  (i  w z ) variation  in  free  space, 
where  K = 1,  such  that,  using  (5), 

u2  + w2  - 1 (13) 

For  the  slot  as  a radiator,  the  wave  in  free  space  must  be  outgoing,  or  moving 
in  the  positive  z direction,  which  requires  that  w be  the  positive  root  of  (13). 
Using  (3),  we  also  have  in  free  space 


ex  (u,  z)  • w h (u,  z) 


(14) 
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Therefore,  since  y (u) 


is  defined  as  h(u,o)/ex(u,o) , 


y(u) 


i / /TT? 


(15) 

Since  we  require  the  positive  root,  we  have  i times  the  root  of  the  absolute 
value  when  u2  is  greater  than  unity,  which  implies  that  y(u)  then  lies  on 
the  negative  imaginary  axis.  Thus  the  imaginary  part  of  the  free-space 
aperture  admittance,  Ya0.  is  negative. 


In  the  presence  of  the  plasma,  (14)  applies  only  in  the  space  above  the 
plasma.  It  will  be  seen,  however,  that  the  plasma  can  be  represented  in  terms 
of  a square  matrix  Q such  that 


e ■ 

V 

X 

= a 

h 

L J 

owter 

h 

(16) 


where  ex  implies  ex  (u,  z)  and  similarly  for  h and  where  "outer"  and  "inner" 
refer  to  either  boundary  of  the  plasma,  the  inner  boundary  being  at  the  ground 
plane  where  z = o.  In  view  of  (14)  applying  outside  the  plasma  and  of  the 
continuity  of  ex  and  h,  the  left-hand  side  vector  is  an  arbitrary  constant 
times  the  vector  [l'].  The  admittance  y(u)  is  the  ratio  of  the  components  of 
the  right-hand  side  vector,  so  solving  (16)  we  have 


y(u)  ~ Can  ~ ~ 3/a  ) 


(17) 


It  should  be  obvious  by  now  that  y (u)  is  an  even  function  of  u,  so  the 
integral  in  (12)  can  be  written  as  twice  the  integral  for  only  positive  values 
of  u.  Also,  it  is  convenient  to  change  variables  in  the  integral  by  letting 
x * ua/2. 


* °°  ^ 

^ - — J y(zx/&)  (s  <■»  */x ) Jx 


(18) 


If  a is  not  too  small,  a further  simplification  can  be  made  by  ignoring  the 
variation  of  y within  each  lobe  of  the  (sin  x/x)2  function. 


% ~ ir  [ I^si  * 

+ Z /((**+')*/*)  f [~£t]*<Jk] 


(19) 

We  have  numerically  evaluated  the  Integrals  in  this  sum  and  found  them  to  fit 
the  following: 

~ yCC*n+l'; TTVfr  ) 

4^  0.1° yC°)  + £r,  7 r*- fr+f/k -o.ts/Cri-z*)]* 

(20) 


-,-V.  . • Vi  «•„«..  , . 


1 


-19- 


Turning  to  the  program  listing,  (20) is  found  in  lines  630,  650,  790,  840,  and 
880.  The  series  is  truncated  by  the  operator  according  to  the  discussion  in 
section  3.1.  However,  the  effect  of  the  omitted  terms  can  be  estimated  ana- 
lytically for  Ya  if  n is  sufficiently  large  at  the  point  of  truncation, 
and  this  correction  factor  is  given  at  line  980  of  the  program.  Eq.  (17) 
is  found  at  lines  640,  810,  and  870.  (Note  that  the  Q matrix  is  represented 
in  the  program  by  the  linear  array  Q.) 

Hie  noise  power  generated  by  the  plasma  and  crossing  the  aperture  is 
the  sum  of  the  contributions  of  all  plasma  layers.  Assuming  Kirchoff's  law 
applies,  the  contribution  of  a given  plasma  layer  is  the  total  power  absorbed 
by  that  layer  when  the  power  flowing  out  from  the  aperture  is  equal  to  that 
which  a black  body  would  emit  if  it  had  the  temperature  of  that  plasma  layer 
and  the  area  of  the  aperture.  The  power  emitted  by  a black  body  at  temperature 
T in  the  frequency  band  df  is,  according  to  Planck's  law, 

p __  A k Ta  df 

nr  (c  T*/‘r-  , ) (2i) 

where  A is  the  area  in  our  system  of  nondimensionalized  units,  k is  Boltzmann's 
constant  and  TQ  is  proportional  to  the  frequency.  At  10  GHz,  To  is  0.48  °K;  but 
we  are  interested  in  thousands  of  °K  temperatures  in  the  plasma.  Therefore 
the  exponential  function  is  quite  accurately  expressed  by  its  linear  approximation 
and 

P - A k T df  / (2 IT  ) 


For  reasons  similar  to  that  which  led  to  this  linear  function  of  temperature, 
the  noise  power  in  an  electronic  circuit  is  equal  to  k Tjj  df,  where  TN  is  the 
noise  temperature.  Using  this  fact  to  remove  k df  from  the  equations,  we 
now  have  the  noise  temperature  at  the  aperture  equal  to  the  sum  over  all 
plasma  layers  of  the  noise  temperature  of  the  layer  as  seen  at  the  aperture. 
Also,  this  layer  noise  temperature  is  the  power  absorbed  by  the  layer  when 
the  aperture  output  power  is  A/ (2 V) times  the  plasma  temperature.  But  the 
linearity  of  the  plasma  absorption  implies  that  the  aperture  output  power 
may  be  set  equal  to  a constant  if  the  power  absorbed  by  the  layer  is  scaled 
by  AT/ (2  Tr'  ) divided  by  that  constant.  In  our  approximation  of  the  aperture 
field  distribution  the  aperture  power  is  ^AjEopGg  where  Ga  is  the  aperture 
conductance.  Thus  the  aperture  noise  temperature  is 


A P 

v\Be\ 


where  A?  is  the  power  absorbed  in  the  layer  when  the  aperture  radiates  and 


has  electric  field  E0  • 


The  power  absorbed  in  a plasma  layer  is  equal  to  the  drop  from  the  inner 
to  the  outer  boundary  in  the  power  flow  P.  Using  the  Fourier  transform 
representations  for  E and  H, 

f * ± i>  £«.  £ fff  e (“j  *)  V,z)  e*r  [ 1 r'J  ' 'M Ju  ^ v dxj 
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where  b ia  the  length  of  the  slot,  (even  though  the  model  is  stated  In  terms  of 
a semi-infinite  slot  having  a uniform  field,  the  noise  temperature  is  proportional 
to  its  length.)  Assuming  that  the  x integration  may  be  performed  first,  its 
result  vanishes  unless  u = v.  Basic  mathematics  requires  that  the  resulting 
function  weights  the  integral  over  all  v as  follows: 


f f (v)  ] exp  [i(u-Ox]dx  dv/  = Z.TT  -f  (u) 

-to  -to 

Therefore 

«*  ^ 

P s TT  t J Re.  h (u,z)\du 

-o* 


(25) 


(26) 


Now  use  (16)  to  give  ex  and  h,  except  that  Q represents  propagation  from  the  inner 


plasma  boundary  to  the  one  of  current  interest  rather  than  the  outermost  one. 
Also, at  the  inner  boundary  use  h (u,  o)  * e0  y 

is  given  by  (11),  and  y represents  y (u) , which  is  given  by  (17) 


, where  eQ  represents  e (u,o),  which 


(27) 


The  result  of  this  equation  is  the  following  form: 

= *TrG  [-%•  j~  (sinx/x)*  [T,  +■  f 

0 l 

+ T3  Ir»(y)  + | rl*  ]dxj  (28) 

where  y is  y (2  x/a).  The  part  of  this  expression  enclosed  by  { } is  similar 
to  (18),  so  it  is  actually  a by-product  of  the  aperture  admittance  calculation. 

We  use  (20)  to  accumulate  the  integral  numerically.  In  the  program  the  array  TA 
represents  T , ...,  T^  and  the  integral  is  accumulated  at  lines  680  and  910. 

The  factor  outside  the  brackets  is  included  at  tine  1010  after  the  iterations  are 
concluded.  Of  course  a is  CAYA.  However  the  aperture  length  b i3  taken  as  27T 
times  the  body  radius,  CYL. 

Far  field  effects  are  calculated  between  lines  1030  and  1200.  The  theory 

for  this  calculation  is  based  on  the  Fourier  integral.  The  form  of  the  fields 

outside  the  plasma  in  terms  of  their  transforms  e and  h may  be  gotten  by  letting 

K«1  and  substituting  the  Fourier  integrals  in  (3)  and  (4).  Also,  in  (3)  we 

recall  that  the  transforms  of  the  fields  vary  as  exp  (iwz).  Thus  in  free  spac«, 

e « - uex/w  and  we  have 
z 

Ex(x/2)^  / (29) 


E (*  z)  * - J («/*) 

-*  (30) 

where  z - d is  the  boundary  of  the  outer  edge  of  the  plasma  and  where  w is  the 
positive  root  of  (13).  We  require  the  limits  of  the  expressions  when  the  distance 
z approaches  infinity.  In  these  limits  the  integrals  have  the  form 
**  to 

[ f(u)  t.  *<Ja  , where  the  phase  t tends  to  infinity. 

Not  only  does  the  phase  tend  to  infinity  but  also  its  rate  of  variation  with  u 
tends  to  infinity  for  most  values  of  u.  The  contributions  to  the  Integral  vanish 
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in  all  regimes  of  u where  the  phase  variations  are  infinite,  and  the  only  contri- 
bution is  from  the  neighborhood  where  does  not  vary.  Thus  integrals  of 
this  type  are  evaluated  by  the  method  of  "stationary  phase".  Now 


<p  =»  ux  4-  )/l  - u3'  (z-d) 

- y - U (z  - d )/<J  ! - U * 
d u 

Therefore  the  phase  is  stationary  when  u *=  uQ,  where 

U = * /r  = Son  e 

O 


(31) 

(32) 


(33) 


where  r is  the  radial  distance  from  (o,d)  to  (x,  z)  and  & is  the  incidence  angle, 
or  the  angle  between  the  line  of  sight  to  (x,  z)  and  the  normal.  The  Taylor 
series  expansion  of  ^about  uQ  is 


0-P-(«-Uo)*'r/(2CoS*e)  + ••• 


(34) 


Use  only  these  terms  in  the  integral  and  replace  f (u)  by  f (uq) , giving 
/ f (m  ) e ‘ * 'f  ( *»)  e'*"  J"exp  c.oS*e')]du 

- oo  6 

= f caS  9 'fJrr^TTy 

(35) 


Therefore* 


t r 

Cx(sir'&Jd)e  cos  &■  JlTT/C  it) 

E2  (xj  z)~  - £x  ( «yz)  ~ta.r>  & 


(36) 

(37) 


Thus  the  far  field  is  transverse  and  is  given  by 

E ~ (sln£>)  d)  e ' r JiTr/ir  ’ (3g) 

Now  (16)  gives  ex  at  z ■ d in  terms  of  the  Q matrix  and  the  values  of  ex  and  h 
at  z = o;  then  (17)  gives  h in  terms  of  the  0 matrix  and  ex  since  y (u)  is  h/ex 
at  z = o;  and  finally  ex  (sin©,  o)  is  given  by  (11).  Therefore 

,‘r  f <3a.x  ~ (3,2  Q,,  *1  Ep  Sin(  j-a  sin  9) 

LQx*  ~ J .Site 

(39) 

where  the  Q matrix  is  evaluated  for  u = sin  ©.  This  is  the  final  equation  for  the 
far  field  of  the  slot  in  the  presence  of  the  layered  plasma  and  in  terms  of  the 
assumed  field  EQ  in  the  slot. 

* 

Throughout  this  report  the  symbol  means  "asymptotically  equals",  usually 
in  the  far  field  limit. 
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Of  principal  interest  in  the  problem  of  plasma  effects  is  the  far  field  at  a 
given  angle  9 relative  to  what  its  value  would, be  in  the  absence  of  the  plasma. 
Thus  the  total  signal  attenuation  is  determined  by  the  ratio  of  the  results  of 
(39)  wi ch/wi thou t plasma.  It  will  be  seen  that  E depends  on  whether  or  not 
plasma  is  present  and  that  the  Q matrix  is  the  identity  in  the  absence  of  plasma. 
Therefore 


Q/,13  **  ~ $ i a ^2  1 

£a(re) 

- Qu 

c 

(40) 


where  superscript  zero  denotes  values  in  the  absence  of  plasma.  In  this  equation 
only  the  factor  involving  the  Q matrix  is  a function  of  the  angle  0,  so  this 
factor  represents  the  far  field  pattern  distortion  function,  which  is  calculated 
in  terms  of  dB  loss  in  lint?  1140  through  1170  of  the  program.  Of  possible  additional 
interest  is  the  absolute  power  density  in  the  far  field.  But  since  the  power 
at  the  aperture  relative  to  the  primary  input  power  is  a function  of  the  circuit 
behavior  of  the  antenna,  which  is  determined  elsewhere,  consider  the  power  gain 
relative  to  the  aperture  power.  This  "aperture  gain"  function  is  the  ratio  of 
Tt\E(rj  /^0/a< ij- 


& = 


a. 

3 

^11  " & '3  ^ a / 

* ^ 

j.  s i * e 

0-2  3 ~ & A /COS  & 

(41) 


The  quantity  is  calculated  in  dB  at  lines  1100  through  1180  of  the  program.  The 
part  of  it  which  does  not  depend  on  the  Q matrix  is  also  output  at  line  1190. 

Note  that  this  latter  auantity  only  approximately  represents  the  "aperture  gain" 
in  free  space,  it  being  necessary  to  correct  it  by  the  ratio  of  the  aperture 
conductance  with/without  plasma  in  dB.  However,  it  is  convenient  to  output  it 
in  this  form  and  a simple  calculation  based  on  the  printout  of  aperture  admittance 
gives  the  correction  factor. 

3. 2. 1.2  Subroutine  PROP 


The  purpose  of  this  subroutine  is  to  accumulate  the  values  of  the  propagation 
matrix  Q and  the  four  temperature  factors  appearing  in  (28).  In  the  calling 
sequence,  N Is  the  number  of  plasma  layers,  DIEL  is  the  array  of  their  complex 
dielectric  constants,  X is  the  array  of  their  thicknesses,  M is  an  input  index 
which  bypasses  the  temperature  factors  if  M is  negative,  V corresponds  to  the 
input  value  of  u or  sin©  , and  Q is  a linear  complex  array  constituting  Q,,  , 

Q,4  , Qa/,  and  QAx  in  that  order.  Error  indication  is  made  by  returning  V equal 
to  the  negative  of  the  value  which  was  input. 

The  equations  from  which  the  propagation  matrix  is  derived  are  given  above. 
In  particular,  substitution  of  the  Fourier  integrals  in  (3)  gives 


dh 

dz 


1 Ke. 


(42) 
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Also,  consider  a given  plasma  layer  In  which  the  complex  dielectric  constant  K 
is  constant.  Then  the  Fourier  Integrals  substituted  In  (5)  give 


(43) 


where  a * 

y s k - u 


Of  course  the  solution  of  (43)  is 

, i rz  -l  rz. 

h * e +■  a*  e 

where  a^  and  a2  are  arbitrary  constants.  Now  (42)  gives 


(44) 


(45) 


c 

X 


i 


- A,  € 


yz 


(46) 


Of  course  the  origin  of  z is  arbitrary  In  all  of  these  equations  so  let  us  place 
It  at  the  inner  boundary  of  the  layer  In  question.  Then,  letting  z = o in  (45) 
and  (46),  solving  for  aj  and  an<*  substituting  back  in  (45)  and  (46),  we 
get  Che  matrix  equation 


N 

■ ^ y 

where 

fcosC  Yz) 

i C r/k)jcn(*-2) 

P = 

1 

(k /Y)Scn(Yz) 

CeS( Yz) 

(47) 


(48) 


Of  course  z in -this  equation  is  the  thickness  of  the  layer.  The  meaning  of  "outer" 
and  "Inner"  in  (47)  refers  to  the  individual  plasma  layer  while  (16),  which  defines 
Q,  refers  to  the  entire  plasma.  Now  both  e and  h are  required  to  be  continuous 
at  all  boundaries  between  layers,  so  the  Q matrix  for  the  first  two  layers  would 
be  the  product  ?2  where  P,  denotes  the  P matrix  for  the  j th  layer  from  the 
Innermost  plasma  boundary,  in  general  the  Q matrix  is  the  cumulative  product 
of  the  P matrices  for  ail  layers,  where  the  addition  of  a layer  In  outward  order 
gives  Q as  P times  the  old  value  of  Q. 

In  terms  of  the  program  listing,  the  value  of  Q is  initialized  as  the  identity 
in  lines  1360  - 1390.  The  effect  of  all  plasma  layers  is  accumulated  in  the  loop 


4 


I 
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fromline  1400  to-  line  H70."  CA  Is  Y at  line  1410.  The  distinct  elements 
of  P are  calculated  asTC8*,CB  and  CF  in  lines  1420  - 1560.  The  cumulative  matrix 
multiplication  PQ  is  done  at  lines  1570  - 1620.  Note  that  the  error  guards 
at  lines  1430  and  1730  - 1760  are  required  by  the  limited  exponent  (10^  38)  0f  the 
L66  computer  word. 

Comparison  of  (26),  (27)  and  (28)  shows  that  Tj  through  T^  are  each  equal  to 
the  layer  temperature  times  the  drop  across  the  layer  in  the  value  of  the  correspond- 
ing combination  of  the  elements  of  the  cumulative  Q matrix  as  given  in  (27).  The 
arrays  TB  and  TC  denote  these  four  combinations  of  elements  of  Q,  TB  being  for 
the  inner  edge  of  the  layer  and  TC  for  the  outer.  As  1340  and  1350  indicate  all 
combinations  are  initially  zero  except  for  the  second,  which  is  initially  unity. 
Current  values  of  TC  are  calculated  at  lines  1640  - 1690.  As  noted  above,  the 
array  TA  is  the  four  T's  in  (28),  and  their  values  are  accumulated  at  line  1710. 


3. 2. 1.3  Subroutine  CRITS 

In  long  experience  with  calculations  of  em  wave  propagation  in  plasma,  we 
have  found  that  spurious  effects  can  be  induced  by  approximating  a smooth  plasma 
distribution  in  terms  of  a series  of  uniform  layers.  Such  spurious  effects  are 
greatest  when  one  or  more  of  the  layers  is  at  critical  electron  density.  They 
tend  to  be  minimum  and  of  no  effect  on  the  results  when  the  critical  density 
is  midway  between  the  densities  of  adjacent  layers.  In  other  words  if  the 
critical  density  point  of  the  smooth  profile  lies  at  a boundary  point  between 
adjacent  layers  of  the  layered  approximation,  the  latter  tends  to  be  a good 
approximation.  This  subroutine  is  designed  to  produce  this  effect,  together 
with  the  main  program. 

In  the  calling  seauence  the  array  FI  is  the  input  profile  of  relative 
electron  density,  XI  is  the  input  profile  distances,  N is  the  number  of  points 
in  the  input  profile,  and  the  array  A represents  NA  (up  to  two)  values  of 
relative  electron  density  which  are  to  be  considered  critical.  (Although  not 
done,  for  convenience,  in  this  application,  it  is  possible  to  consider  the  point 
where  7*  vanishes  in  addition  to  where  K does  as  a critical  point.)  During 
execution  NA  is  changed  to  indicate  the  number  of  locations  used  in  the  array 
XA.  The  meaning  of  the  array  XA  is  a set  of  points  in  the  coordinates  of  the 
input  profile  at  which  layer  boundaries  must  be  placed.  Layers  of  constant 
thickness  may  be  used  between  all  values  of  XA. 

The  subroutine  is  designed  to  work  under  the  assumption  that  the  input 
plasma  profile  has  only  one  relative  maximum.  If  it  has  more  than  one,  sub- 
seouent  crossings  of  the  critical  points  will  be  ignored.  The  logic  of  the 
routine  can  be  seen  by  studying  the  listing. 


3.2.2  SLOC 

Fig.  7 is  a listing  of  this  program,  which  consists  of  a main  program  and 
five  subroutines.  Subroutine  CRITS  is  identical  to  the  same  subroutine  in 
SLOP,  so  it  is  not  described  in  this  section. 


3. 2. 2. I Main  Program 

The  first  58  lines  of  this  program  are  similar  to  the  first  55  lines  of 
SLOP,  where  the  arrays  of  coi^lex  dielectric  constants  and  thicknesses  of  the 
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plasma  layers  are  defined.  The  only  exceptions  are  the  lines  240,  360  and  370. 
The  array  YIP  is  the  ratio  of  the  thickness  of  the  layer  to  the  cylindrical 
radius  at  its  Inner  edge  and  CAYR  becomes  the  cylindrical  radius  of  the  outer 
edge  of  the  plasma. 


The  aperture  admittance  and  plasma  noise  temperature  are  calculated  in 
lines  600  through  1230.  The  derivations  are  based  on  the  assumption  of  a 
cylinder  of  radius  R having  a gap  of  width  a about  its  circumference.  The 
plasma  is  composed  of  a series  of  co-axial  layers,  each  of  uniform  but  differ- 
ing density.  In  cylindrical  coordinates  let  E^m  0,  Ezm  0 if  /*|>a/2,  and 
E = E if  Jzj^ta/2,  on  the  cylinder  of  radius  R.  Thus  throughout  all  space 
tfie  magnetic  field  H = , Ep  = 0 and  all  fields  are  independent  of  0.  In 

the  cylindrical  coordinate  system  0 Z ) Maxwell's  equations  are 


'B  £> 
B z 

U± 

B z 


*2 


(49) 

(50) 

(51) 


Where  K is 
separation 
exp  (i  u z 
transforms 


a function  of  f3  only.  These  equations  have  the 
of  variables  and  superposition  as  in  the  planar 
) separates  instead  of  the  same  function  of  x. 
are 

0°,  Z)  * _S  Cry  u)  exp'iu*)J* 


(»'  • 


/ i' 


I 

OZ 


same  solutions  by 
case  except  that 
Typical  Fourier 


(52) 

(53) 


and  similary  for  E^and  H. 

Let  us  preserve  the  convention  implied  by  the  derivation  used  in  the 
planar  case,  that  the  aperture  admittance  is  positive  for  power  flowing  outward 
through  the  aperture.  Then  Y is  defined  as  the  negative  of  H/Eq  in  the 
aperture  and  the  analog  of  (9^  is 


- J Ez  («  z> 

-oo 


(54) 


Substituting  the  transform  for  H,  reversing  the  order  of  integration,  and  using 
the  inverse  transform  of  Er  we  get 

*- 

i IfJ*  X = -*rr  f h (r/ «)  zz  Cr  u)<Ju  (55) 

-oe 

Now  ez  (R,u)  is  identical  to  ex  (u,  o)  as  given  by  (ll).  Therefore  (12)  and, 
as  a result,  (18),  (19)  and  finally  (20)  apply  for  the  cylindrical  as  well  as 
the  planar  case.  Of  course,  however,  the  admittance  function  y(u)  la  the  neg- 
ative of  the  magnetic  field  divided  by  the  tangential  electric  field  at  radius 
R for  a cylindrical  wave  which  varies  as  exp  (i  u z) . 
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Substitution  of  the  Fourier  transforms  gives  Maxwell's  equation  as 
follows,  on  elimination  of  the  f>  component. 


r*h  - cK 

(56) 

(57) 

where  X Is  given  by  (44).  Elimination  of  h gives 

37  ( 77s-)  + K r = 0 

(58) 

In  free  space,  where  K = 1 and  w,  this  is  Bessel's  equation  of  zero  order 
in  the  variable  pw . The  solution  corresponding  to  outgoing  radiation  in  our 
convention  is  HQ  (w^>),  which  is  the  Hankel  function  of  the  first  kind.  Also, 
(56)  states  that  h is  i/w  times  the  derivative  of  this  function.  Therefore 
for  free  space. 


y (u) 


i Ht  (wR) 
w H0  (wR) 


(59) 


where  we  recall  that  w ■ J 1 - u^  . Even  when  plasma  is  present  the  external 
variations  of  ez  end  h are  the  same,  so  let  R represent  the  cylindrical  radius 
of  the  outer  edge  of  the  plasma  and  use  the  matrix  Q to  represent  the  plasma 
effect  on  the  tangential  fields.  Therefore  let 

ft  iv  jl 

Al  J = \ 

(60) 

where  A is  an  arbitrary  scalar.  Now,  since  y ( u ) ■ - h/ez  at  R, 


w 4**  H0("*r)  + i 4,1 


(61) 


Turning  to  the  interpretation  of  the  computer  program  listing,  note  that 
in  calls  of  subroutine  HANK  the  variables  A,  B,  C and  D are  the  real  and 
imaginary  parts  first  for  Hq  and  second  for  H . Therefore  (59)  shows  up  in 
lines  680  and  910,  where  CA  and  CB  are  w H and  1 H^,  respectively.  Similarly 
(61)  is  seen  in  lines  720  and  970.  The  lo§ic  also  includes  a branch  covering 
imaginary  values  of  w,  which  uses  the  identities 


t hr,  Cx) 

H/(L  *)  * - C*/n)  k,(n) 


(62) 

(63) 
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where  K and  K ^ are  the  modified  Bessel  functions.  Noting  that  w equals 
i times°a  real  quantity,  called  W in  this  branch  of  the  program,  and  that  calls 
of  HANKI  give  K as  A and  as  B,  (59)  shows  up  in  line  1020.  Also,  (61) 
is  represented  8y  line  1090. 

The  derivation  of  the  plasma  noise  temperature  is  the  same  as  before  up 
through  (23).  The  power  flowing  out  across  a given  layer  boundary  at  radius 
f>  is 

OC 

p=  - 7 T f>  Re{  JXf  &z(Ku)  CKf>[c(  U-V)ZJ 


® ^ (f;  « ) h*Cr;U)]ju  (64) 

Now  the  fact  that  ez  and  h are  transformed  by  the  matrix  Q from  their  values 
at  f = R,  where  h = - y ez  and  ez  = eQ,  gives 

Rt  [ ez  U )]  ~ I > 

— Rz(Q,t< 0^*  + Q i/  0*  ) Re  ( y) 

+ ~ Q„  ) I^(y)  4-  ) Irl*] 

(65) 


Therefore 


\ { i f (sin  x/k)*  [t  + rx  i~Cr) 

^ O 

- 7}  R'  Cx?  + t¥  lylx]  dxj 


(66) 


As  before  the  array  TA  represents  Ti,  ...»  T^  in  the  program.  The  integral 
is  accumulated  at  lines  760  and  1130;  and  the  multiplying  factor  is  included  at 
line  1230,  where  a is  CA/A  and  R is  CYL. 


The  far  field  effects  are  calculated  between  lines  1250  and  1470.  In  the 
far  field  ez  is  the  solution  of  Bessel's  equation  for  n ■ o and  for  argument 
w /=>  , so 

z)  ■=  if  IV  A,  c • <X  *)du 


f ^V"  + u 

00  " 


where  A is  a function  of  u determined  by  the  boundary  conditions  and  where  the 
asymptotic  form  of  the  Hankel  function  has  been  used.  Now  the  phase  is 
s tationary  In  this  integral  when 


a » W,  = 2/ r 


= SOI  8 


(68) 
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where,  as  before,  9 is  the  angle  between  the  line  of  sight  and  the  normal. 

(The  fact  that  now  r Is  taken  from  the  cylinder  axis  Instead  of  from  the  edge 
of  the  plasma  makes  negligible  difference  in  the  limit  of  infinite  r.)  The 
phase  expands  as  In  (34)  about  this  point  of  stationary  value  and  therefore  (35) 
still  applies.  The  result  is 

~ * A cos  6 1 V/*  (69> 

where  A Is  evaluated  for  u = sin  ©.  Now  (49)  and  (50)  give,  at  this  point  of 
stationary  phase 

Ep(rj9)~  -td/nO  (70) 

Therefore  the  far  field  Is  transverse  and  is  given  by 


E (r^o)  ZA  t L*‘//'r 

(71) 


The  coefficient  A is  determined  by  (60)  from  elimination  of  h at  p ■ 
E (r  a)  ~ (u-Z/z)  CQ„  Q*t.  ~ 


TT  « IT  (i  Q^H0  - qu  H,) 


R.  Thus 


(72) 


where,  of  course,  u ■ sin® 
functions  of  argument  wR  . 

R = R,  so  p 

P 


, w » cos  & , and  Hq  and  Hi  are  the  Hankel 
In  the  absence  of  plasma  Q Is  the  Identity  matrix  and 


£ (*■_,&) 

X (w  1 

e0 

E'C^o) 

- ^2*  +■  <-  Qu  H, 

e.° 

(73) 

The  part  of  this  equation  within  [ J is  at  line  1390  In  the  program  and  line 
1400  gives  Its  value  In  dB  down.  The  "aperture  gain"  Is  477  r2  times  the  square 
of  the  far  field  divided  by  the  aperture  power. 


c = ** 

T-  • _ 2 

1 st/ 1 y\ 

Q"  Q 2*  ~ Q 13.  Qa! 

RGj.  (tt*)3- 

L “*/*.  j 

+ L &'Shi t/\*/ 

(74) 

The  value  of  this  quantity  In  dB  Is  the  FORTRAN  variable  C at  line  1450.  Except 
for  the  lack  of  normalization  by  the  free-space  conductance,  as  noted  in  SLOP,  the 
FORTRAN  variable  D Is  the  free-space  value. 


3. 2. 2. 2 Subroutine  WFP 

The  purpose  of  this  subroutine  is  to  accumulate  the  values  of  the  prop- 
agation matrix  Q and  the  four  temperature  factors  appearing  In  (66).  In  the 
calling  sequence,  LAYERS  Is  the  number  of  plasma  layers,  DIEL  Is  the  array  of 


* 


. 


v.  V:  ■ 


-37- 


t 
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their  complex  dielectric  constants,  DK  is  the  array  of  their  thicknesses, 

YIP  is  the  array  of  their  relative  thicknesses  in  terms  of  the  thickness 
divided  by  the  inner  radius,  CAYR  is  the  radius  of  the  innermost  boundary,  U 
is  the  input  value  of  u or  sin  & , and  Q is  the  linear  complex  array  represent- 
ing Qn,  Q^2  Q21  and  Q22*  Error  indication  is  made  by  returning  U equal  to  the 

negative  of  its  input  value. 

In  a given  plasma  layer,  having  a constant  value  of  K,  (58)  becomes 
Bessel's  equation  of  zero  order. 


t'z/Ctr)  + e2  = o (75) 

Where  * denotes  differentiation  with  respect  to  Also  (56)  may  be  written  as 

k * */r)  (76) 

Hopefully  not  confusing  the  reader,  let  us  denote  a pair  of  linearly  independ- 
ent solutions  of  (75)  by  u and  v.  Then  in  the  given  layer 


€z  =•  a,  u -h  V (77) 

A = (i  k/r)  Cs.,  v ')  (78) 

Now  let  subscript  o denote  the  value  of  a quantity  at  the  inner  boundary  of 
the  layer.  Then  the  undetermined  coefficients  a and  <12  can  be  solved  for  in 
terms  of  eZQ,  h0,  u0,  v0)  u'  and  v/  . Now, in  order  to  simplify  the  calculations 
and  since  we  use  nonstandard  functions  (see  subroutine  PROP)  for  u and  v, 
let  v0  = ii0'  » o.  Then  the  vector  (ee,  h)  is  the  propagation  matrix  P times 
(ezo»  ho),  where 


~ c * */C  * *a' ) 

Lku'/Cruj  v'/< 

Subroutine  PROP  outputs  u/u0>  v/v  ‘ , u7u0,  and  v1/  v^as  P,QQ,  R,  and  S 
at  line  1770.  CA  represents  T at°line  1710. r Lines  1800  through  1850  are  the 
accumulation  of  the  Q matrix  from  multiplication  by  P. 

Lines  1860  through  1950  and  line  2010  involve  the  accumulation  of 
the  four  temperature  factors  in  (66).  Inspection  of  (65)  in  comparison  with  that 
equation  reveals  the  reason  for  the  form  of  each  term  in  lines  1880  - 1910. 

These  factors  are  corrected  by  the  ratio  ^/R  in  lines  1920  and  1940  because 
(64)  contains  while  (66)  has  R.  Finally,  the  difference  is  taken  in  terms 
of  "outer"  minus  "inner"  (TC  - TB  in  line  1950)  to  give  the  power  lost  in  the 
layer,  because  of  the  negative  sign  in  (64). 


3. 2. 2. 3 Subroutine  HANK 


This  routine  calculates  the  real  and  Imaginary  parts  of  the  Hankel 
functions  of  order  zero  and  unity  and  for  real  argument.  In  the  calling 
sequence,  X is  the  argument,  FJO  is  the  real  part  of  Hq,  FYO  is  the  imaginary 
part  of  H , FJ  is  the  real  part  of  H^,  and  FY  is  the  imaginary  part  of  H^. 

The  method  uses  the  polynomial  approximations  on  pages  369-370  of  Abramowitz 
and  Stegun  (Ref.  1).  These  polynomial  equations  are  written  in  the  form 
of  continued  products  in  order  to  obviate  the  calculation  of  powers. 


3. 2. 2. 4 Subroutine  HANKI 

This  routine  calculates  the  modified  Bessel  function  of  the  second 
kind  of  orders  zero  and  unity.  In  the  calling  sequence,  X is  the  argument, 
FKO  is  K , and  FK  is  K..  The  method  is  the  polynomial  approximations  on 
pages  378-379  of  Ref.  I.  Again  the  continued  products  form  is  used  for 
economical  computation  of  power  series. 
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3.2.2. 5 Subroutine  PROP 

The  purpose  of  this  routine  is  to  calculate  the  elements  of  the  matrix 
P ao  in  (79).  It  is  more  general  than  needed  by  subroutine  WFP  since  the  vari- 
able EN  in  the  calling  seauence  refers  to  the  order  of  Bessel's  equation.  The 
variables  XZR  and  XZI  are  the  real  and  imaginary  parts  of  the  value  of  Yp 
for  the  inner  boundary  of  the  layer,  and  Y is  the  relative  layer  thiokness. 

The  outputs  are  C,  D,  CP,  and  DP,  which  represents  u/uQ,  v/v^  , u//\iQ,  and 
v;/  v0 ' respectively. 

The  method  uses  series  expansions  in  the  relative  layer  thickness  Y,  where 
Y is  the  thickness  divided  by  the  radius  at  the  inner  boundary.  The  fact  that 
u/uQ  and  being  on  the  diagonal  of  P,  must  approach  unity  as  Y approaches 

zero  implies  that  the  leading  term  in  series  expansions  must  be  unity.  Similarly 
the  lowest  order  term  of  the  other  two  quantities  is  at  least  first  order. 
Therefore  let 

OO 

* / + x>^  Z.  Cj  (80) 

v r / 

+■  jZj  ) (81) 

oa 

u'/ua  — O + O C ; 

° J--,  J (82) 

OO 

v'/v/  - / + ZL  (j  + i )■=(.• 

° j*i  J (83) 

In  these  equations  X is  the  value  of  *7®  at  the  inner  boundary  of  the  layer 
and  the  Cj  and  d.  are  each  eoual  to  a complex  constant,  which  depends  on  j, 
times  the  power  of  Y.  The  Cj  and  di  in  (82)  and  (83)  are  identical  with 

those  in  (80)  and  (81)  since  the  second  two  quantities  are  the  derivatives  of 
the  first  two  with  respect  to  XY.  Now  u/uQ  and  v/v'  must  satisfy 


+ * TJ  +(z*-r,J)«~o  (84) 

where  z»>7»in  the  layer.  In  particular 

z = x Ci  + y)  <85> 

where,  of  course,  X is  constant  and  Y is  considered  variable.  Substitution 
of  (80)  and  (81)  in  (84)  and  separation  of  terms  of  like  powers  of  Y,  keeping 
in  mind  the  implicit  proportionality  between  cj  or  dj  and  Y^ , gives  recurrance 
relations  for  the  cj  and  d,.  The  same  relation  holds  for  both,  so  let  aj 
denote  either  cj  or  dj 

J(j  + ')<Ej  + [(J-O*- niJyA 
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This  equation  applies  liberally  for  either  term  as  long  as 


>j  “ ° * 


d_j  = o 


1/ (XY) 


It  is  not  difficult  to  see  these  equations  in  the  listing  of  the  computer 
program.  However,  see  ref.  2 for  a more  detailed  account  of  the  theory, 
programming,  and  checkout  of  this  routine,  which  is  identical  in  that  report 
with  the  present  routine.  Reference  2 also  contains  additional  details 
concerning  subroutine  WFP. 


3.2.3  SLOS 

Fig.  8 is  a listing  of  this  program,  which  consists  of  a main  program 
and  five  subroutines.  Subroutine  CRITS  is  identical  to  the  same  subroutine 
in  SLOP,  so  it  is  not  described  in  this  section. 


3. 2. 3.1  Main  Program 

With  two  exceptions  the  first  60  lines  of  this  program  are  similar 
to  the  first  58  lines  of  SLOC.  One  exception  is  that  the  aperture  con- 
ductance GA  is  not  initially  set  to  zero  since  the  logic  is  not  a function 
of  this  quantity  as  an  input.  The  other  exception  is  that  the  array 
DIEL,  which  is  defined  at  line  590,  is  the  square  root  of  the  complex 
dielectric  constant  of  the  layer. 


The  aperture  admittance  and  plasma  noise  temperature  are  calculated 
in  lines  610  through  1120.  The  derivations  are  based  on  the  assumption 
of  a sphere  of  radius  R having  a gap  of  width  a at  a given  latitude.  The 
plasma  is  composed  of  a series  of  concentric  spherical  layers,  each  of 
uniform  but  differing  density.  In  spherical  coordinates  (r,  © , <p)  let 
°*  E#“  0 if  > a/  (2R) , and  E*  - E0  if  /© \<  a/ (2R)  on 

tne  sphere  of  radius  R,  where  SQ  is  the  location  of  the  center  of  the 
aperture.  Thus  throughout  all  space  H • H#  , E*  ■ 0 and  all  fields  are 
Independent  of  <Z> . Maxwell's  equations  are 


CtE:e)  ~ ■=  L r H 

^ [h  s ih  e)  = - t Mr  r £Tt  si*  & 

hOM)  = 
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* Figure  8 (continued) 
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where  K is  a function  of  r only.  Elimination  of  the  components  of  E from 
these  equations  gives 

(sir®  sime))  + Kr&(£??(*H))+K**H  = 0 ^ 

Separation  of  variables  can  be  accomplished  by  letting 

'ieCjke  = + 

(95) 

which  is  Legendre's  equation  of  order  unity  and  degree  n.  Thus  let 

HCr,*-)  - 2.  Pn'(coS&) 

h - i iyb7 

Now  (93)  implies  that  E^  has  this  same  dependence  on  &. 

PC 

Ea(rjQ)  ^ 2.  %Crjn)  p/fcos  B)  (97) 

h*t 

Also,  (92)  implies  that  Er  depends  on  Pn  (cos  & ). 

~ -21  er  n)  Pn^cos9)  (98) 

h - ^ 

To  calculate  the  aperture  admittance  start  with 

7 r 

A I? , l*  y&  = IT  Ea(R;s)  sc*  * de 


where  A is  the  area  of  the  slot. 

, cl*  Jtr 

A = a ir  ft  J Jt'*  & do 

~ jht 

» 

“ *Tr*K[cos(*.-j %-)  +£-)] 

= airR^(x2  -X,)  (100) 

Substitute  (96)  and  (97)  in  (99)  and  evaluate  the  Integral  using  the 
orthogonality  of  the  Legendre  functions  with  respect  to  the  degree  n. 
Then,  using  y (n)  ■ h (R,  n)/  e&  (R,n) , 

, jr  h{t)+  l)  I gfl  /*  , , 

V = (z/t+O  l~£Tl 


**.  -X, 

Where  eQ  ■ etf  (R,  n).  Orthogonality  gives 


e CXFWt)<ix 

0 + j r^cxjdx 


(101) 


(102) 
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where  x^  and  ^ are  implied  by  (100) , 


SubstiCuCion  of  the  series  solutions  in  (93)  and  (94)  gives 

1 K r v W-  (r  h) 

K r ^[-TT  HF"  (r  h)]  + “ n (n+1)J  h * 0 


(103) 


(104) 


If  K is  constant  this  last  is  the  equation  of  the  spherical  Bessel  functions 
of  order  n and  argument  Kr.  In  free  space  the  solution  corresponding  to 
outgoing  waves  is  the  spherical  Bessel  function  of  the  third  kind,  hn  (r) . 
(From  here  on  we  drop  ' ' from  the  notation  since  our  time  convention  makes 
it  clear  that  we  need  only  those  functions  which  tend  to  e *r/r  at  large  r.) 
Thus  y (n)  in  free  space  is  given  by 


y(* ) -= 


~ ln/R 


(105) 


In  the  presence  of  plasma  let 


= a 


(106) 


Replacing  R by  R_,  the  radius  of  the  outer  edge  of  the  plasma,  (105)  gives 
the  ratio  h/etfatP that  radius.  Then  (106)  gives 


/0)=  - 


where 


Z = c»/*f  -t 


aoy'i 


(108) 


Now  in  terms  of  the  FORTRAN,  XX  at  line  640  is  (x£  - ) and  EON(M)  becomes 

eQ  /E0  for  M ■ n after  line  800.  At  line  920,  B is  the  coefficient  of  y (n) 
in  the  sum  of  (101).  CC  at  line  930  is  the  denominator  of  (105)  and  line 
940  completes  the  accumulative  value  of  Ya  • CB  at  line  950  is  Z as  given 
by  (108),  and  CA  at  line  960  is  Ry  (n)  as  §iven  by  (107).  Yfl  is  accumulated 
at  line  970. 

• 

The  power  flowing  through  a sphere  of  radius  r is 

? = -ir  r*  f 2L  2L  Kt  **(*)£* 

— t hi  */  A*/ 

***2-  eJ1-1 '■,('■;*>  <109> 

As/  Lin  + ij 


■>  > "*  V>  • 


*•'  ' 


Now  use  (106)  to  get 
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Kef  h*  etfO,  n)  in* (r,  n)]  - r /ej*  [ (4az  Q,+) 

+ Re.  (4s,Q,t  +■  Q3AQ,f)  Re(rty) 

■b  r„  ( 4*  a Q*  ~Q*lQu  ) Tr*  (fry) 

+■  Re  ) ! R J (no) 

where  Ry  Is  Ry  (n)  as  given  by  (107).  Now,  therefore,  (23)  gives 
00 

* ~ % Zr,  7*n  +0  1^1  h + rz  Rs(*y) 

+ Tj  r*(f?y)  +■  tv  (in) 


This  equation  appears  in  lines  980  and  990,  except  that  G0  is  included  at  line 

1120.  a 

The  asymptotic  form  of  hn  (r)  is  - ie  ir/(i  n r).  Therefore  in  the  far 
field,  H varies  as  e lr/r.  Now  (93)  implies  that  Ee«  H in  the  far  field 
and  (92)  implies  that  r2  | Er}2  tends  to  zero  as  r approaches  infinity.  Thus 
the  far  field  is  transverse  and  can  be  related  directly  to  either  E&  or  H. 

£„(>■<,>  ~ - U(Rr>‘''>  r>,  a) 

cr  i”  r"  s) 

n 

(112) 

The  quantity  Z,  which  is  given  by  (108),  is  the  free-space  impedance  (ratio 
of  ^(Rp.  n)  to  h (R  , n))  at  radius  R , so  etf  - Zh  at  the  outer  edge  of 
the  plasma.  Using  (i06)  we  then  get  P 


A (ftp,  n)  - 


g,  (<*„  Qaa  ~ 

z4„  - "r**, 


(113) 


The  "aperture  gain"  is  2 TTr2  times  lE-  (r,  fi>)  I 2 divided  by  the  total 
aperture  power,  fcA  | Ee  | 2Ga.  1 


G-  - 


4 /r  e,CrJe>)/£0  /* 

R*  (**-*,)% 


(114) 


In  the  program  the  quantity  Qi(R  , n)#c  hR  (Rp)  )J  is  stored  versus  n in  the 
•rray  HN  at  line  1010.  The  value  of  this  quantity  in  free  space,  which  is 
gotten  by  replacing  Rp  by  R,  by  letting  Q be  the  identity  matrix,  and  by  re- 
placing Z by  the  reciprocal  of  y (n)  as  given  by  (105),  is  stored  in  HN0  at 
line  1020.  The  sum  in  (112)  for  each  case  is  accumulated  by  neglecting  to 
keep  account  of  the  net  phase,  as  at  lines  1240,  and  1250;  i.e.,  the  effect 


I 
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of  1 "n  is  gotten  by  adding  the  new  term  without  It  to  1 times  the  previously 
accumulated  sum.  This  Is  simply  a FORTRAN  short  cut. 


3. 2.  3. 2 Subroutine  WFP 


The  purpose  of  this  subroutine  Is  to  accumulate  the  values  of  the  propa- 
gation matrix  Q and  the  four  temperature  factors  appearing  in  (111).  In  the 
calling  sequence,  LAYERS  is  the  number  of  plasma  layers,  DIEL  is  the  array  of 
the  square  root  margin  of  their  complex  dielectric  constants,  DK  Is  the 
array  of  their  thickness,  YIP  is  the  array  of  their  relative  thicknesses, 

CAYR  is  the  radius  of  the  innermost  boundary,  N is  the  value  of  the  order  n 
of  spherical  wave,  and  Q is  the  linear  complex  array  representing  Q2^ 

and  Q22.  Error  indication  is  made  by  returning  N - 1 - n.  In  other  words 
N < 1 Indicates  error  since  n > 0. 

Let  prime  denote  differentiation  with  respect  to  a variable  z ■ J K r . 
Then  in  a uniform  medium  (103)  can  be  written  as 


(115) 


Also  (104)  is 

z2  (rh)"  + [z2  “ n (n  +1)J  (rh)  - 0 


(116) 


Let  u and  v denote  two  linearly  independent  solutions  of  (116).  Then  in 
the  layer  in  question 


rh  ■ Sj  u + a2  v (117) 

e#  ■ -1  (a^  u'  + a2  v')  (118) 

As  in  section  3. 2. 2. 2 let  ve  “ uc'  ■ 0.  Then  the  vector  (rh,etf  ) equals  P times 
the  vector  (r0  h0,  e^  ),  where 


<*/u„ 


-dr/z)CK,/u,') 


t z. 


V/V 


(z./l)  ('''/**’) 


(119) 


This  matrix  is  generated  in  lines  1540  - 1630  and  the  accumulation  of  Q by 
multiplication  by  P is  done  in  lines  1640  - 1690.  The  four  temperature 
factors  in  (111)  are  accumulated  in  lines  1700  through  1780  and  1840.  Each 
term  in  Q , including  r,  which  la  called  A at  line  1720,  can  be  seen  at  lines 
1730  - 1760. 
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3. 2. 3. 3 Subroutine  HANK 

This  routine  calculates  the  spherical  Hankel  functions  of  various  orders. 

In  the  calling  sequence,  M is  the  number  of  orders  to  be  calculated,  X is  the 
argument,  and  HN  is  the  complex  array  of  Hankel  functions.  The  function  of 
order  n is  stored  in  the  (n  + l)t*1  location  of  HN.  The  method  uses  downward 
recurrence  starting  with  arbitrary  values  at  a sufficiently  large  order  to 
get  the  real  parts.  Then  the  imaginary  parts  are  generated  by  upward  recurrence 
starting  with  the  closed-form  values  for  orders  zero  and  unity.  The  real  parts 
are  normalized  so  as  to  agree  with  their  two  lowest  order  closed-form  values. 

This  routine  is  similar  in  theory  ana  development  to  subroutine  HANK  in  ref.  2, 
except  that  that  routine  is  for  ordinary  Hankel  functions. 

3. 2. 3. 4 Subroutine  LEG 

This  routine  calculates  the  associated  Legendre  polynomials  of  order  unity 
and  various  degrees  in  the  calling  sequence,  X is  the  argument,  cos  0;Y  is  sin  6;  M Is 
the  number  of  degrees  to  be  calculated;  and  the  array  P is  the  output  Legendre 
functions.  As  in  HANK,  the  function  of  degree  n is  stored  in  the  (n  -I-  1)  th 
location  of  the  array  P. 


pie  method  uses  upward  recurrence  starting  with  P©  (x) 
and  p|(x)  ■ sine  (line  2400).  The  recurrence  relation  is 

r\  ( *)  =■  (z*  + 0x  PnU?  - (n  + 0 P*_/0 


o (line  2390) 


(120) 


which  is  line  2460. 

3. 2. 3. 5  Subroutine  PROP 

This  routine  is  identical  in  function  and  calling  sequence  to  the  routine 
of  the  same  name  in  SLOG,  which  is  described  in  3. 2. 2. 5.  The  listing  differs 
from  the  former  one  only  at  lines  2590,  286C  and  2870.  These  differences  occur 
when  (80)  and  (81)  are  substituted  in  (116)  instead  of  in  (84).  Then,  instead 
of  (86),  the  recurrence  relation  is 

j( +*i (j-oyt-j.,  + - »> fr+o] 

•+  x* (&j-z  + y *- j -3  + - o (121) 

Therefore  n^  is  replaced  by  (n?  + n ),  as  in  line  2590.  The  variable  E at 
line  2860  represents  the  negative  of  the  coefficient  of  aj_i  divided  by  the 
coefficient  of  aj.  The  variable  F at  line  2870  is  the  same  kind  of  ratio 
with  reapect  to  Sj_2‘ 
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Figure  9 Listing  of  Program  ABCD 
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3.2.4  ABCD 

Fig.  9 is  a listing  of  this  program,  the  objective  of  which  is  to  transform 
the  results  of  calculations  for  conditions  at  the  aperture  to  conditions  at 
the  antenna  input.  The  method  assumes  that  the  antenna  can  be  represented  as 
a linear  two-port  network  and  that  sufficient  experimental  data  exist  to  define 
the  network  parameters  which  relate  voltages  and  currents  at  the  input  to 
electric  and  magnetic  fields  at  the  aperture. 

See,  for  example,  chapter  1 of  reference  3 for  a more  detailed  discussion 
of  two-port  networks.  Such  networks  are  describable  in  terms  of  the  following: 


c-h  P 

a + 8 y* 

(122) 

Ap  - 

BC  - 1 

(123) 

Where  Y is  the  admittance  at  the  input;  Ya  is  the  aperture  admittance;  and 
A,  B,  C and  D are  complex  constants  which  do  not  change  when  external  effects 
change  the  values  of  Y and  Ya.  Rather  than  attempt  to  calculate  A,  B,  C and 
D from  some  physical  model  of  the  antenna,  let  us  rely  on  experimental  data 
to  infer  their  values,  assuming  that  our  previous  calculations  of  Ya  are 
correct  First  we  have  the  calculated  aperture  admittance  in  free  space, 

Ya  , for  which  let  the  measured  input  admittance  be  Y^.  This  gives  us  one 
eqSatlon  in  the  four  unknowns  using  (122).  A second  equation  is  already 
given  by  (123).  A third  is  obtained  if  Y,  is  the  input  admittance  which  is 
measured  when  the  aperture  is  covered  tightly  by  a sheet  of  metal  foil.  Such 
shorting  of  the  aperture  implies  Ya  -oesince  the  electric  field  must  vanish. 

Then  (122)  gives  Y2  ” D/B.  One  more  equation  is  needed  since  there  are  four 
unknowns.  To  get  this  last  equation,  consider  the  effect  of  placing  a thin 
resistance  sheet  over  the  aperture.  Regardless  of  the  assumed  geometry  of 
the  aperture,  its  admittance  is  always  calculated  from  a linear  combination  of 
wave  admittances.  The  boundary  conditions  for  a thin  resistance  sheet  are 
that  the  tangential  electric  field  is  the  same  on  either  side  and  the  tangential 
magnetic  field  has  a discontinuity  equal  to  the  current  in  the  sheet.  This 
boundary  condition  implies  that  the  wave  admittance  is  equal  to  its  value  in 
the  absence  of  the  thin  sheet  plus  the  conductance  of  the  sheet.  In  other  words 
the  sheet  is  in  parallel  with  the  aperture.  Then,  to  a good  approximation,  the 
aperture  admittance  in  the  presence  of  the  sheet  is  Ya<j  + Ya,  where  Yg  represents 
the  conductance  of  the  sheet. 

Now  the  three  simultaneous  equations,  which  supplement  (123), are 


y,*  + /,y^B  - c- yiop~ 

Kl  8 - P « 0 

y3*  + y3(%o+ys)t-Z-()i0+K)t>*0 


(124) 

(125) 

(126) 
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Substitute  D from  (125)  in  (124)  and  (126)  and  divide  each  equation  by  B, 
giving  two  simultaneous  equations  in  A/B  and  C/B,  the  solution  of  which  is 


a/s  - (y3  - y*')  ys  s/y  - y, ) -yt, 

•=-/e  = (A/e/y  + (y,- 

Now  divide  both  sides  of  (123)  by  B^,  using  for  D/B,  to  give 


(127) 

(128) 


B - I /sJ(A/B)y^  - C/B 


(129) 


This  value  of  B then  gives  A and  C from  their  ratios  to  B,  which  are  already 
known,  and  D ■ Y2B.  In  the  program  the  names  A and  C are  temporarily  used 
for  A/B  and  C/B  at  lines  150  and  160,  which  represent  (127)  and  (128).  Line 
170  is  (129)  and  the  final  values  of  A,  C and  D are  derived  at  lines  180  - 200. 

If  a calibrated  signal  is  input  on  a transmission  line  of  characteristic 
admittance  Y0  and  the  complex  voltage  reflection  coefficient  R is  measured, 
the  total  voltage  at  the  input  is  (1  + R)  times  the  input  voltage.  Also,  the 
total  current  at  the  input  is  Y0  (1  - R)  times  the  input  voltage  and  therefore 
the  admittance  Y at  the  input  is 


(130) 


This  forms  an  alternative  way  of  prescribing  input  conditions  for  the  above 
calculations  of  A,  B,  C and  D parameters.  Thus  the  complex  voltage  reflection 
coefficients  for  each  of  the  three  conditions  can  be  used  to  specify  Y^  Y2  and 
Y3,  given  Yq.  This  is  done  using  (130)  in  lines  70  - 120  of  the  program. 

Of  course  Inversion  of  (130)  gives 

k - ^ 

yo  * y (131) 


which  is  done  at  line  140  if  the  Y values  were  input.  The  dB  loss  by  re- 
flection in  free  space  is  calculated  at  line  220. 

Now  if,  owing  to  the  presence  of  plasma,  the  aperture  admittance  takes 
on  some  new  value  Ya,  then  (122)  gives  the  new  value  of  input  admittance  Y 
and  (131)  gives  the  reflection  coefficient  at  the  input.  This  gives 


R - — 

K CA  + F yz)  f 


(132) 
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This  equation  is  found  at  lines  260-270,  and  lines  280-300  calculate  the 
amplitude,  phase  and  dB  loss  corresponding  to  R. 

Plasma  effects  on  radiation  patterns  were  calculated  in  SLOP,  SLOC  and  SLOS 
with  respect  to  the  total  electric  field  at  the  antenna  aperture,  E0,  which  is 
itself  affected  by  plasma.  The  matrix  equation  corresponding  to  (122)  is 


= F 


(133) 


where  V and  I are  the  total  voltage  and  current  at  the  input  port,  Eo  and  Ho 

are  the  electric  and  magnetic  fields  at  the  aperture,  and  F is  a scaler  constant  which 

converts  the  aperture  fields  to  equivalent  voltages  and  currents.  Invert  this 

equation,  keeping  (123)  in  mind,  let  V * V0(l+R)  and  I * Y0vo(l”R) » and  use 

(132)  for  R.  Then 


= 


* K K /F 


(134) 


where  VQ  is  the  input  signal  voltage.  Using  this  equation  with/without  plasma. 


y.  (*  + * -hc-+  py* 

K ca  + b £)  + c +-  p yA( 


(135) 


This  quantity  is  calculated  at  line  310  and  in  terms  of  dB  loss  at  line  320. 

Finally  the  noise  power  at  the  aperture  needs  to  be  transformed  to  the  input.  We 
calculate  this  effect  by  continuing  the  physical  model  of  equating  the  emission 
to  the  absorption.  That  is,  if  the  power  input  to  the  antenna  were  set  equal  to 
the  noise  temperature  which  is  seen  at  the  aperture,  as  calculated  above,  then 
the  power  at  the  aperture  would  equal  the  noise  temperature  which  is  seen  at  the 
antenna  input  port.  In  symbols,  let  the  aperture  noise  temperature  Ta“%  -0/V0/ 2 
while  letting  the  noise  temperature  at  the  input  port  rH  \Ff0\ » 


Tn  * 


*3.  *Ta. 

\ y»  ( *+  & + c + Ja 


(136) 


This  equation  is  applied  at  lines  330-340,  where  TN  denotes  Ta  and  T denotes  Tn. 
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4.  Conclusion 

The  models  presented  in  this  report  are  the  lowest-order  examples  of  more 
general  theory  for  plasma-covered  aperature  antennas  in  canonical  geometries. 

The  possible  generalizations  are  of  various  types.  Let  us  discuss  some  of 
these  here. 

The  infinite  integrals  for  the  aperture  admittance  and  noise  temperature  induced 
by  the  plasma  in  the  planar  and  cylindrical  geometries  are  a subject  for  gener- 
alization in  terms  of  numerical  techniques.  The  techniques  used  here  are  best 
adapted  to  the  situation  in  which  the  aperture  is  large.  In  particular  the  inte- 
gral in  (18)  was  evaluated  by  neglecting  the  variation  of  y (2x/a)  within  each 
lobe  of  the  sin  x/x  function;  i.e.,  Y (u)  is  assumed  to  vary  little  when  u 
varies  by  2V/a,  which  can  be  true  only  when  a is  large.  Obviously,  this  numer- 
ical technique  is  easily  improved  upon  by  evaluating  the  integral  over  each  lobe 
of  the  transform  of  the  aperture  field  distribution  by  the  trapezoidal  rule, 
for  example.  With  little  revision  of  the  logic,  one  could  even  let  the  operator 
truncate  these  integrations  as  is  now  done  with  the  sum  over  the  lobes  of 
sin  x/x.  Thus,  in  a given  iteration  of  the  trapezoidal  approximation,  we  would 
add  the  effect  of  including  all  points  midway  between  those  already  included, 
via  a do-loop. 

The  assumed  uniform  field  in  the  aperture  could  be  changed  to  some  other,  fairly 
arbitrary,  distribution.  The  only  restriction  which  suggests  itself  is  that  the 
electric  field  parallel  to  the  edge  of  the  aperture  must  vanish  at  the  edge. 

The  general  techniques  used  in  this  report  would  still  apply.  The  near  field 
effects  (aperture  admittance  and  plasma  noise  temperature)  are  gotten  by  invoking 
conservation  of  energy  at  the  aperture  between  the  assumed  field  distribution  and 
the  separatlon-of-varlables-solutlons  of  Maxwell's  equations  just  outside  the 
aperture.  And  the  far  field  is  gotten  in  relation  to  the  assumed  aperture 
electric  field,  ignoring  matching  to  the  magnetic  field.  An  interesting  example 
of  an  aperture  field  distribution  is  the  case  in  which  Eg  varies  as  cos  0 
in  the  slot  on  a sphere.  Then  the  field  does  not  vanish  at  © = 0 or  7T , 
as  it  does  when  the  aperture  field  is  uniform;  but  in  general,  the  field  depends 
on  both  9 and  0. 

The  present  models  are  for  one-dimensional  apertures  (infinitely  long  for  planar 
geometry  and  circumferential  for  cylindrical  and  spherical  geometries).  It  is 
possible  to  revise  them  to  treat  two-dimensional  apertures.  The  simplest 
example  would  be  a rectangular  aperture  in  which  the  electric  field  is  parallel 
to  the  short  dimension  and  has  a cosine  distribution  of  amplitude  along  the  long 
dimension.  In  the  calculations  of  aperture  admittance  and  noise  temperature, 
the  single  Integral  or  sum  would  be  replaced  by  a double  integral  or  sum.  Far 
field  effects  would  be  calculated  in  only  slightly  different  ways  from  the 
present  models.  For  example  in  the  planar  case,  the  present  model  gives  the 
relative  effect  of  the  plasma  on  the  far  field  for  directions  which  are  in  the 
plane  of  the  short  dimension  of  the  aperture.  All  that  needs  to  be  added  is 
the  effect  of  plasma  on  the  far  field  in  directions  for  which  the  magnetic  field 
lies  in  the  plane  of  propagation.  Circular  apertures  in  a ground  plane  can  be 
analyzed  also.  The  relative  far  field  is  the  same  as  for  a rectangular  aperture 
in  terms  of  the  polarization  of  the  antenna  relative  to  the  plane  of  the  line 
of  sight.  The  near  field  effects  for  a circular  aperture  require  reformulation 
in  terms  of  the  transforms  of  the  aperture  field  distributions. 

It  is  possible  to  generalize  these  models  to  Include  more  than  one  antenna 
aperture.  The  theory  is  quite  similar  to  the  above,  except  that  the  field  distri- 
bution at  the  body  surface  is  that  for  the  array  Instead  of  for  a single  aperture. 
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The  transform  of  this  field  distribution  is  slightly  more  complicated,  and  the 
admittance  is  a symmetric  matrix  rather  than  a scalar.  If  the  same  plasma 
profile  covers  all  elements  of  the  array,  the  admittance  matrix  comes  from  a 
straightforward  calculation  and  far  field  effects  for  each  aperture  are  no 
different  than  they  would  be  for  a single  aperture.  Additional  pattern  distortion 
may  be  induced  by  the  plasma  effect  on  the  admittance  matrix,  and  such  effects 
could  be  calculated  using  the  concept  of  an  n-port  network,  for  example. 

An  array  has  the  capability  to  produce  a narrow  beam,  which  will  be  distorted 
by  the  plasma  even  if  the  plasma  is  uniform  over  all  elements  of  the  array. 

This  effect  comes  about  because  the  far  field  loss  induced  by  the  plasma  is  a 
strong  function  of  angle  of  incidence.  Thus  the  parts  of  the  beam  nearer  the 
normal  are  less  attenuated  by  plasma  than  are  those  farther  from  the  normal, 
which  tends  to  pull  the  beam  toward  the  normal.  This  effect  can  be  calculated 
either  by  modeling  the  beam  shape  and  applying  the  attenuation  by  plasma  at 
a number  of  closely  spaced  angles  in  the  beam  or  by  analytically  relating  the 
beam  pointing  error  to  the  rate  of  change  of  attenuation  with  incidence  angle. 

The  case  in  which  the  plasma  varies  from  one  aperture  to  another  can  not  be 
calculated  directly  by  these  canonical  models  since  their  essential  feature 
is  the  assumption  that  the  plasma  varies  only  in  the  normal  direction.  However, 
one  can  estimate  the  effects  of  such  transverse  gradients  by  treating  each 
aperture  as  a separate  canonical  problem.  The  far  field  effects  can  be  derived 
by  considering  the  effect  of  differing  plasmas  on  the  coherent  addition  of 
all  elements  in  forming  the  beam.  Such  a calculation  can  be  done  by  analytical 
approximation  in  terms  of  the  rate  of  variation  of  the  effects  of  plasma  on 
amplitude  and  phase  versus  distance  across  the  afray.  The  near  field  effects 
can  be  approximated  by  letting  the  mutual  admittance  between  two  given  elements 
be  equal  to  the  average  value  obtained  when  both  are  covered  by  either  plasma. 

All  of  these  effects  can  be  calculated  for  the  case  when  the  plasma  has  random 
fluctuations  among  array  elements  as  well  as  when  the  variations  are  smooth. 

Some  of  these  model  improvements  already  exist  in  another,  much  older,  computer 
program  which  has  never  been  documented  formally.  However,  the  logic  and 
numerical  techniques  in  that  program  are  cumbersome,  so  we  intend  to  incorporate 
such  improvements,  when  needed,  in  the  present  programs. 
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